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Abstract 



The thermal and density corrections, in terms of the isospin chemical potential /ij, to the mass of the pions, 
the decay constant and different condensates are studied in the framework of the SU (2) low energy effective 
chiral lagrangian at finite temperature in the two phases: The first phase \fii\ < m and the second phase 
|/ii| > TO, being m the tree-level pion mass. As a function of temperature for /.jj = 0, the mass remains quite 
stable, starting to grow for very high values of T, confirming previous results. However, there are interesting 
corrections to the mass when both effects (temperature and chemical potential) are simultaneously present. 
At zero temperature the tt^ should condense when ji^ = ^m^. At finite T, the condensed pion acquires a 
thermal mass in such a way that a mixture, like in a superfluid, of a condensed and normal phase appears. 
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Preface 



There is a reasonable agreement that Quantum Chromodynamies (QCD) is the right theory that explains 
strong interactions. The QCD degrees of freedom are quark fields and gluon gauge fields, which possess color 
charge and. up till now, have not been detected as isolated free particles. The theory cannot explain the 
confinement of the quarks and gluons inside the hadrons, and it is not possible to obtain, in general, analytical 
solutions, except for processes at high energies or "hard processes" (also equivalent to small distances) where, 



due to the asymptotic freedom, perturbation theory is possible |GW73H |GW73a| , |GW74 |Pol73| . Non- 
perturbative methods like the Operator Product Expansion (OPE) and QCD Sum Rules [$VZ79t or the 



Nambu Jona-Lasinio (NJL) effective model [NJL6f | include corrections to the perturbative analysis, with 
quark-gluons degrees of freedom, in terms of condensates, responsible for the hadronic resonances at low 



energy. Numerical simulations of QCD on the Lattice | tH78 , tH81|, on the other side, turn out to be a real 
non-perturbative approach to the theory, starting from first principles. 



Chiral Effective field theories, like the Linear Sigma Model [GML6C], the Non-Linear Sigma Model | Wei6S 



and Chiral Perturbation Theory (xPT) [Wei79, GL84, GL85], have been a strong and useful tool for under- 
standing QCD at low energies or "soft processes" , since they involve hadronic degrees of freedom preserving 
the symmetries of QCD and allow us to establish a connection with the concepts of Current Algebra (the 
precursor of QCD). Among all hadrons, pions play a special role in the dynamics of hadronic matter because 
they are the lightest hadrons and therefore are produced at a higher rate in high-energy reactions. The 
discussion of pion dynamics is a precious window to explore the non-perturbative structure of QCD. 

In the last 20 years, there has been a special interest in hot and (or) dense processes in scenarios like 
the early universe, relativistic heavy ion collisions, and inside neutron stars. High energy processes in a 
thermal bath have been studied in different frames, searching for new phenomena and states of matter like 



vn 



Chiral Symmetry Restoration (see for example | NRZ96{ ), the Quark-Gluon Plasma |Mul85, KS04, Kod04 



and confinement of quarks and gluons into liadrons, Color Glass Condensate [ILMOl, FILM02, McL04|, etc. 
In all these cases we deal with different phase transitions. A main problem in this context is to find the 
appropriate order parameters, in terms of observable quantities associated to the behavior of a quasi-particle 
in a medium. 

The dependence of the pion properties on finite temperature has been studied in a variety of frameworks, 



such as thermal QCD-Sum Rules |DFL9C], Chiral Perturbation Theory at one loop | GL87b and two loops 
[lSch93| , |Tou97|] , the Linear Sigma Model [|Lar86i |CL90i |DLR94| , the Mean Field Approximation |lPis82 



BCDC^92| , the Virial Expansion [ LS90 , Sch91 |, etc. There seems to be a reasonable agreement that the pion 
mass m-^{T) is essentially independent of T, except possibly near the critical temperature Tc, where m-^{T) 
increases with T, and that the pion decay constant ,f-n{T) vanishes at the critical temperature. 

The introduction of in-medium processes via isospin chemical potential has been studied at zero temper- 
ature in xPT pSOl| , |KT01| , ^TV02a| in both phases ^ m^) at tree level. 

In-medium properties at finite density have been discussed for a variety of phenomena, as for example. 



the chiral condensates |PLL_03, PLLW03 |, the anomalous decays of pions and etas | CRK03 |, etc. As the 
density increases, both the quark condensates and the decay rates diminish. 

Interesting results concerning the structure of the QCD phase diagram, including temperature effects, 
have been achieved for the case when we have simultaneously baryon chemical potential and isospin chemical 
potential. Two completely different approaches have confirmed a qualitative change in the phase diagram 



as soon as the isospin chemical potential starts to grow [TK03, BPRC03|. These problems concerning the 
structure of the phase transition diagram, have also been handled in the frame of two color QCD in four 
dimensions (a QCD-like model) | |STV02l| as well as in three dimensions [DN03]. The problem with baryonic 



chemical potential has been considered in the frame of xPT [AEGN95| and also using the finite pion number 



chemical potential AAA02 |. 

Different properties of QCD (or QCD-inspired models) under these kind of circumstances have also been 
analyzed in the lattice approach. For QCD with three and two colors, extensive work has been carried out 
[iKTSOll , |KTS02| , |KS02b| , |KS024 |KTS03|| in connection with the behavior of different order parameters in 
several phase transitions in the /i/T plane as, for example, the transition to a diquark phase, the chiral 
condensate, etc. 



Isospin asymmetry docs exist in nature in the case of neutron stars. Since 1971, the idea of a pion 
condensate in the core of neutron stars has been considered in connection with the coohng process of a 



neutron star (see for example |YP04]). This idea is a motivation to study the behavior of pions at extreme 
isospin densities, searching for phase transitions and considering the possibihty of having such scenario, 
isospin asymmetric regions, in RHIC and ALICE experiments. EventuaUy, these considerations are relevant 
in early stages of the universe. 

The main results of this thesis concerns the discussion of thermal radiative corrections, including finite 
isospin chemical potential role on both pion phases, enlarging in this way the existing literature on the isospin 



chemical potential effects in the pion phase transitions |STV02b, KS04]. The results of this thesis with respect 



to the thermal masses and decay constants in the first phase can be foimd in [LV03] and the behavior of the 



thermal masses in the second phase in [ LV04 1 . The non-tivial evolution of pion masses and decay constants 
may have also important phcnomcnological consequences for the diagnosis of the pion gas. in the central 
rapidity region of relativistic heavy ion collisions, when looking for signals as dilepton or photo-production. 

This thesis is organized as follows: In chapter |l|, I will introduce Chiral Perturbation theory, the renormal- 
ization procedure involved, and the finite temperature and density formalism. In chapter ^ the derivation 
of the effective Lagrangian, the propagators necessary for radiative corrections, and the expansion criteria to 
solve the problem of the non-diagonal propagators in the second phase will be presented. In chapter ^, the 
self-energy corrections, decay constants, and condensates in the first phase will be derived. Chapter ^ and |^ 
will be devoted to the calculation of the self-energy corrections and condensates in the regions > m and 
/ii 3> m. Finally, in chapter ||, I will summarize the results and conclusions, presenting also an outlook. 
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Chapter 1 



Introduction 



This chapter is a brief introduction about the theories and tools that will be used to compute radiative 
corrections within a thermal bath for the case of a pure pion gas; i.e., xPT and Thermo-Field Dynamics. 



1.1 Chiral Perturbation Theory 

xPT is a successful effective field theory that makes use of the chiral symmetry, which is an intrinsic symmetry 
of massless fermions. Since light quarks (u, d) possess masses of a few MeV (the s quark can also be handled 
in this frame), it is reasonable to think of them as massless, allowing a separation of the quark fields in right- 
and left-handed components. Then, it is possible to construct an effective Lagrangian with hadronic degrees 
of freedom that preserves chiral symmetry. The idea of the construction of an effective theory in this context 



was first proposed by Weinberg |Wci79 



The chiral symmetry breaking is usually parametrized by a non- vanishing pion decay constant /tt. using 



Partial Conservation of Axial Current (PCAC) |GML60, Nam60|. According to the Goldstone theorem 



[GSW62|, a consequence of the spontaneous symmetry breaking will be the appearance of massless bosons. 
In the case of xPT, where the SU{N)ji x SU{N)]^ chiral symmetry grou reduces to a SU{N)y, N"^ — 1 
massless bosons appear, according to the number of generators of the group SU{N) that do not annihilate 
the vacuum. In the case of SU{2), for the two lightest flavors u and d, the massless goldstone bosons are 
associated with the three pion pseudo-scalar fields. In fact pions are pseudo-goldstone bosons, due to the 



intrinsic mass of the quarks. 

In the same way as it happens in the theory of superconductivity, we expect that quarks should form 
Cooper pairs of bound quark-antiquark states, due to strong attractive interactions, speciaUy if the quarks are 
masslcss. The vacuum state with a quark pair condensate is characterized by the order parameter {qq) ^ 0, 
which also breaks the chiral symmetry, allowing the quarks to acquire effective masses as they move through 
the vacuum. Then the pions acquire a mass due to the explicit symmetry breaking induced by the small 
quark masses and also due to the existence of this quark condensate. 

One of the advantages of xPT, for example in the case of two flavors, is that we can forget about 
the isosinglet scalar field a because it is heavier than the pseudo-scalar fields [nia ~ (400 — 1200)A/ey in 
comparison with m^± w XAQMeV). Then, it is possible to concentrate on pion interactions only, in a scenario 
like a pion gas, especially when dealing with finite-temperature corrections. In effective models, like the linear 
cr-model, it is possible to discuss the properties of the cr-meson, and its influence in the low-energy dynamics 



[CL90, DLR94]. Here, I will only concentrate on the pion dynamics, dominant when finite temperature and 
density effects are taken into account. 



A new comprehensive review of xPT can be found in 1 3ch02 1 . 



1.1.1 Chiral Lagrangian 

Let us proceed in the frame of the SU (2) chiral perturbation theory. The most general chiral invariant 
expression for a QCD-extended Lagrangian under the presence of external hermitian-matrix auxiliary fields, 
has the form 

Cqcd{s,p, V, a) = Cqcd + ^qcd + 97^(iV + Tsa^)? - ^(•'^ - il5P)q, (1-1) 

where Cq^jy is the usual QCD Lagrangian with zero mass, Ca is the anomalous part^ and s, p, v, and a are 
scalar, pseudo-scalar vector and axial external fields defined as 

v^, = v'^{x)t''/2, s = sO(x)-f s"(a;)T", 

= a"(a;)rV2, p = p° (x) + p'' ix)T'' 



^In principle, it includes also the anomalies of the fcrmion determinant but it will not be considered here, because it does 
not contribute to the corrections involved in this work 
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where s°, a° arc real functions and t° arc the Pauh matrices. Consider the independent local 

SU X SU (2) L transformations of the left- and right-handed quarks in flavor space (chiral transformations) 



(l + 75)'7->3i?,(l + 75)g, (l - Ts)*? ^ - 75)9, (1-3) 

which induce gauge transformations of the external fields 

v'f, + a'^, = gR{Vf, + af,)g\j + igRdf,gl^, (1.4) 
= 9l{v^- a^,)g{ + igLdf,g\, (1.5) 
s' + ip' = 9R{s + ip)gl- (1.6) 

Then, it is possible to construct a generalized chiral-invariant Lagrangian containing only pion degrees of 



freedom, where the external fields transform in the same way as in Eq.(1.6). The effective low-energy 



Lagrangian will be expressed as an expansion in powers of momentum P on a certain scale A^^, 

C^=C2+C^ + Ce + --- + C^, (1.7) 

where £„ ^ 0{P"') and is the Wess-Zumino-Witten Lagrangian (the anomaly contribution). 

According to Weinberg power counting [ Wei79|| , d,a,v,s,p ^ 0{P). This means that the external mo- 



mentum must be smaller than a certain scale A^ in order to proceed with an expansion in a series of powers 
of momentum^ . 

We will start with the 0{P^) chiral Lagrangian 

£2 = ^Tr[iD^UyD^^U + U^x + X^U] (1.8) 

with 

Df^U = d^U -i[v^,U]~i{a^,U}, (1.9) 

X = 2B{.s + ip), (1.10) 

U = U"^{e''''^^'^/f)U''\ (1.11) 



We will see afterwards that = 47r/ ~ 1.2GeV, where / is the pion deeay eonstant in the ehiral limit 



At this point, / and B arc arbitrary constants with dimension of mass, and U is the vacuum expectation 
value of the unitary-matrix field U . 



The most general ©(P^) chiral Lagrangian has the form 

£4 = \h{Tr[iD^UYD^U]f 

+ ll2Tr [{D^U)^D,U] Tr [{D^U)^D''U] 

+ ^{h + k) (Tr [xU^ + Ux^]Y 

+ \UTr [{D^U)^D^^U] Tr [xU^ + Ux^] 

+kTr [{L^M + UR^,) {U^Lf^" + R^-'U^)] 

+leTr [iL^^Df^UiD^Uy + iR^,{D^U)^ D^U] 

-^h {Tr [xU^ ^Ux^]f 

+ l{h,+h^)Tr[x^^] 



-\{h,-h^)\Tr 



X +r 



Tr [xf - Tr 



-2h2Tr [L^^L^^" + R^.R^'"] , 



(1.12) 



with 



Rf^iy = d^r^ - d^r^ + i[r^, r^], 



(1.13) 
(1.14) 



The li are the original Gasser and Leutwyler coupling constants [ GL84 , GL85 and the hi are couplings 
to pure external fields, and their values depend on the model. For calculations at the one-loop level, it is 
enough to keep terms up to O(P^), as we will see soon. 



We can see that our extended QCD-Lagrangian contains the external fields s, p, v, a, coupled to the 
scalar, pseudo-scalar, vector and axial-vector currents, respectively. Now, it is easy to derive these currents 
and other quantities in the low-energy theory. For example, let us consider the axial current. Using the 
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high-energy QCD Lagrangian, it is given by 



\Qlt^lbT°'q= —SQCD[s,p,v,a], (1.15) 

where S'qcd[s,p, a] is the QCD action. In the low-energy description of QCD, we can obtain this current 
following the same procedure: 

5 S 

where the index n in denote 0{P'"'). 

If we want to calculate condensates, for example the chiral condensate, (also known as quark condensate) 
{qq), the procedure is the same: 

(0|g-g|0) = (0|J,|0) = -{Q\-LsQCD[s.p,v,am. (1.17) 

The chiral condensate in the low-energy regime can be identified with 

5 S 

{qq) = -(0|^|0) = + im)(3) + ■■■ + (1-18) 

and it is proportional to B in the low energy effective theory. 



1.1.2 Power counting, B, f and 

It is obvious that the effective Lagrangian is difficult to handle because it is expressed in terms of an infinite 
series in powers of the fields involved; for such reason we need to truncate this series. Consider the case of 
masslcss QCD without external fields. Proceeding as we did in the previous section, the Lagrangian, the 
axial current and the chiral condensate are (see App. ^) 

£2(0,0,0,0) = ^Tr[dUdU^, (1.19) 

A^i)^(0, 0,0,0) = l-Tr[id^U{^T%U^}-i{^T%U}d,U^], (1.20) 
(gg)(i)(0,0,0,0) = -{0\^BfTr[U + U^]\0). (1.21) 
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Note that this Lagrangian is independent of the vacuum state value. The lower contribution to these quantities 
(tree- level), expanding in pion fields, are 



= (1-23) 

(1.0) 



(ot)(i,o) = -2i?/^ (1.24) 



respectively, where the sub-indices (n, m) mean ©(^"Tr™). 
According to the PCAC relation.'^ 

(0|A;:(x)k''(p)) =^p^/^5''''e-^^'^ (1.25) 
where /^r is the pion-decay constant. In the case of j^^^^, 

(0|A^i,i)^(x)k^(p)) = -!dM'^'^{x)\'K\v)) = ip,f5'''e-^^\ (1.26) 



From eqs. (1.24) and ( 1.26| ) we conclude that / is associated to and Bp is related to the quark 



condensate. In fact / and —2Bp are the pion decay constant and the quark condensate, respectively, in the 
chiral limit. ^ 

Consider now the case of QCD with massive quarks without external fields: Cqc'd{M, 0,0,0), where M 
is the quark mass-matrix (in the SU(2) case, M — diag(jnu,rnd))- The free effective Lagrangian in terms of 
the U fields is 

^—Tr[dUdU^ + 2BM{U + U^)\. (1.27) 



The effective potential is then 



Veff = -^Tr[2BMiU + W)]. (1.28) 



By construction, the vacuum expectation value of the fields U is defined with the fields tt" = 0. A general 



^Stricllj' speaking, the concept PCAC is related to the divergence of the axial current, but actually is known by the relation 



in eq. < \l.2(\ ). 

*If we take into account higher terms in the expansion, it can be proved that their contribution is negligible in the chiral 
limit. 
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expression for U is 



U = cos if + iT ■ n sin tp, 



(1.29) 



where n is a unitary vector and ip, for now, is an arbitrary angle to be derived. It is easy to see that the 
value of U that minimizes the effective potential is U = 1. 

Now, expanding the Lagrangian in terms of the pion fields, we get the free Lagrangian: 



(1.30) 



The constant (m„ + mj,)B is then associated to the pion mass mj. We will refer to the tree-level pion mass 
as = (m„ + mii)B. 



If we compare with eqq. (1.24) and (1.26), we get the well-known Gell-Mann Oakes Renner (G-MOR) 



relation at the tree level |GM0R6J 



■i(TO„ + md){qq) = mlfl 



(1.31) 



This relation opened a new way to our understanding of strong hadron dynamics and chiral symmetry, since 
it relates the the pion masses and decay constants with parameters of the fundamental theory: the quark 
masses and the quark condensate. 



The next terms of the Lagrangian are powers of the form /^(tt^//^)™. If we impose that the vacuum 
sectation value of (tt^//^) ~ P^/A^, wf 
dimensional regularization (see App. |B. j| ) 



expectation value of (tt^//^) ~ p^/A?,, we can neglect higher terms in the Lagrangian. In fact, if we use 



1 

p 



7r"(a;)7r''(a;)|0) 

df^k iS, 



P 



(27r)4 fc2 - m2 + it 



'Jab 



In 



A2 



(1.32) 



where A is the MS divergent term and A is a scale factor associated to the dimensional regularization 
procedure. Forgetting for a moment the divergent term (keeping in mind, however, that the theory is 
renormalizable in an effective sense), certainly the term Inm^/A^ is in the range of 0{P^) for a reasonable 
value of A. We can set then the chiral-scale factor as A^^, = 47r/. Since /^r ^ 93MeV, the theory is valid up 
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to Ax - l.2GeV. 

Now wc can truncate the series according to the number of loop corrections. For example in the calculation 
of 1-loop self-energy corrections, we need to consider corrections up to 0{'k^) in the fields. If they include 
the £2,4 which is of ©(P®), then we must take into account the next contributions of the £2^ up to O(P^) 

£tree — ^2,0 + £2,1 + ^2,2 , (1.33) 

4 

^1-loop = £2,™ + £4.0 + £4,1 + £4,2, (1-34) 

6 4 

C-2-loop = + £4,m + -^e.O + -^6.1 + ^6,2, (1.35) 

m— m—0 

etc. 



I would like to remark at this moment that, if we go to higher orders in the expansion, we will loose 
predictive power. In fact, we need experimental values of some observables to fix the different coupling 
constants, at each order in the expansion. For higher orders we have more and more couplings, since the 
theory is only renormalizable in an effective sense. This point will be discussed more in detail in the next 
sub-section. 



1.1.3 Renormalization 



The chiral Lagrangian is not renormalizable in the usual sense, since it involves an infinite number of couplings. 
Nevertheless, by construction, the different coupling constants involve the infinities that appear in dimensional 
regularization: 

(1.36) 

where the different h and 7^ are tabulated in App. |B.l| 

The renormalization procedure for self-energy corrections works as follows: The corrected propagator is 

Dc = D iDY,D+ ^ D{l + ii:D)-\ 
D^^ = D-^+iS, (1.37) 
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327r2 



T 1 "1" 

li+ln—-X 
A^ 



where D is the free propagator matrix defined as 

D-^^^^ fd\e^P^^-^^^{x), (1.38) 

and S is the self-energy matrix. UsuaUy D and S are diagonal (or anti-diagonal); but we will see that this 
will not be the case in the second phase where \fii\ > m. 



Following the usual renormalization procedure, we re-scalc the fields tt* = ^/Ziir^^ with Zi = 1 + ■ S 



and are of the order of the perturbative parameter P^/A?, = 5^. Since [S] = 2 in mass units and is 



dimensionless we see that S ~ 0{P^S^) and ^ 0{6^). The renormalized propagator is 

= Z,,iD-\p)^,-^{p),,+0{PHl) 

= tD-\p),,-J:n{p)^, (1.39) 



with Zij = ^yZiZj. The value of Z is chosen in such a way that the corrected propagator does not have 
corrections proportional to p^, or, in other words, the coefficient that multiply the term in the corrected 
propagator must be 1. 



As an example, let us consider a free propagator and the self energy correction of the form 

iD-Hp) = (po-fo)'-p'-a', (1.40) 
I](p) = S^ + SVo + Sy. (1.41) 



This expression for S, as we will see, is valid in the first phase. 



Choosing Z = 1 -|- E^, the renormalized self-energy will be 



^ ^"n + ^RPo- (1.42) 
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Expansion in terms of mass corrections 



Unfortunately, in general, the self-energy S will not have the form of eq.(1.41), but could be a complicated 
function of the external momenta, as indeed will be the case for the second phase in the region of high isospin 
chemical potential values. As we want to compute mass corrections, it is possible to expand the self-energy 
in terms of these corrections. In the rest frame, where p = 0, the energy will he po = rrifi = mt + (5m, where 
rriR is the renormalized mass, rrit is the tree level mass and 5m is the correction due to the self-energy terms. 
Then 

I](mfl) = S(mt) + S]'(mt)(5m + il]"(mt)(M^ + ... 

= YP[mt] + Y}[mt]mR + Y,''[mt]m\ + 0{5mf, (1.43) 

where the argument inside the brackets indicates the mass around which the self-energy expansion was 
computed. 

The masses are defined as the poles of the determinant of the propagator matrix at zero 3-momentum, 
i.e. they will correspond to the solutions of 



l^"'(p)|p=o = 0- (1-44) 



As we said before, if the self-energy has a complicated form, we expand the renormalized mass in the rest 
frame in powers of corrections 5m to the tree-level mass. Since the renormalized masses will be extracted as 



solutions of 



D-R^mR) 



0, we only need to compute the 5m corrections. 



PCAC and condensates 

In the case of the PCAC relation, where the axial vector current is saturated with one pion (0|A°|7r^), the 
physical pion is the renormalized pion tt^. The relation will then be 



(0|^;i(7r)|<) = miWZirM). (1.45) 
As the tree-level part of the axial-vector current is proportional to a single pion, the square root of the 
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rcnormalization constant ^fZTa can be expanded in terms of &Zai which is of order by^: 



(1.46) 



The different Za values are set by the self-energy correction conditions, explained before. 

For the case of condensates, the renormalization procedure is not necessary, since all the divergences 
coming from radiative corrections cancel with the li constants. 

1.2 Finite temperature and chemical potential 

The properties of a system in thermal equilibrium can be determined from the grand partition function 



where H is the Hamiltonian of the theory. /3 = 1/T is the inverse of the temperature^, being the Lagrange 
multiplier of the energy, and /i^ are the chemical potentials where Ppn are the Lagrange multipliers of the 
different conserved quantities. With this, the thermal average of an observable is 



In field theory, if we consider the fields 0^, representing the different particles of the theory, the partition 
function is then given as a functional integral 





(1.47) 



(0) = -Tr[e-/^ 



o\. 



(1.48) 




(1.49) 



integrated along an imaginary time path | Mat55 | with 




(1.50) 



^In our unit system, the Boltzmann constant fcs is taken as ks = 1 
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— oo 



oo 







— oo — i3 



-i0 



Fi gure 1.1; Integration path in TFD. The fields that "live" in the imaginary time xq — ^/?/2 are called Thermal Ghosts. 



satisfying the Kubo-Martin-Schwinger (KMS) boundary conditions |Kub57, MS59 



(I>{xq — i/3, x) = (j>{xo,x) for bosons, 
(xo — i(3,x) — —(l){xo,x) for fcrmions. 



(1.51) 



1.2.1 Thermo-Field Dynamics 

To construct the usual generating functional, we can set xq — —t and follow a path of integration in the 
complex time-plane as in Figure 



1.1 



If we take the limit r ^ oo, the fields will be split into two regions; the "real" fields that live on the real 



axis of time, and the so called thermal ghosts. This means that the Fock-space is duplicated | UMT82 |. This 
formalism is named Thcrmo-Ficld Dynamics (TFD). The formal algebraic formulation of TFD was proposed 



by [Oji81|. The path integral formulation for the fermion fields, and the real-time Feynman rules for gauge 



theories with fermions, at finite temperature and density, were developed by [ KSW85 | and NS84 , MOU84 |. 

One of the consequences of the appearance of thermal ghosts is that we will get four two-point propagators 
instead of one for each field, i.e., if we denote by 4>i = 4>{x) the real field and by 02 — (j>{xo — {(3/2, p) the 
ghost fields, we have 

DT,^{x,y-T) = ((O|r,^,(x)0,(y)|O)), (1.52) 
12 



where by |0)) we understand the thermal vacuum state or "populated" vacuum state, and the propagator 



with the real fields Dthh is the so-called Dolan-Jackiw propagator (DJp) |DJ74| 
The prescriptions for loop calculations are: 



• The thermal insertions appear only through loop calculations. 



• Thermal ghosts must be included in calculations at two or more loops 



Since in this work we will deal only with one-loop corrections, we need only to consider the DJp. I will not 



discuss the construction of the DJp, which isbased on the KMS boundary conditions indicated in eq, 1.51 and 



the time-order of the fields along the path described in Fig. 1.1. A general formula for the DJp for bosonic 
fields is 



Dd.j{p) 



- — iim 



D{kQ + irj.p) - D{ko - ii],p) 



nB{po)[D{po + ie,p) - D{po - ie,p)], 



(1.53) 



where —iD is the free propagator in momentum space defined in equation ( 1.38 ) and ns is the Bose-Einstein 
distribution defined as 



1 



pPx _ I 



(1.54) 



For a free propagator of the form 



iD-\p) = - 2apo ~ b, 



(1.55) 



the corresponding DJp, following eq.1.53 is 



Ddj{p) 



p^ — 2apo — b + ie 



27rnB(|po|)(5(p^ - 2apo - b). 



(1.56) 



To distinguish the propagator at zero temperature and the pure thermal part of the DJp, we will refer to 
them as the free-propagator and the thermal insertion, respectively. 



13 



1.2.2 Covariant formalism and finite chemical potential 



Consider the case of a ferniionic field within a dense medium, asymmetric with respect to baryon number. If 
the baryon number density is n = ■0^-0, the introduction of a chemical potential into the Lagrangian density 
as a Lagrange multiplier of this conserved quantity gives 



£ + /in = 'ijj{i"fd — m)ip + fi^j''^tp ~ ■0(17(9 + 70/i — to)0 
= 0(i7^[5^ - i/i(5^o] - m)0. 



(1.57) 



This means that the chemical potential acts like a gauge field [ Wel82 , Act85 1 . The corresponding propagator 
is 



^(7P + 7oM + "^) 



(1- 



The presence of a thermal bath breaks Lorentz invariance. since it plays the role of a privileged reference 
frame. Nevertheless, formally we can "restore" the Lorentz symmetry in such a way that it is possible to give 
a covariant form to the equations, through the introduction of the fourt-velocity between the observer 



and the thermal bath [ Wcl82 1 . In the rest frame (i.e. where the thermal bath is at rest with respect the 
observer), we have 

u= (1,0,0,0). (1.59) 



Although this work is performed in the rest frame, we use the symbol 
is possible to use a covariant formalism. 



3^0 



recalling that in principle it 



1.2.3 QCD at finite isospin chemical potential 



As was indicated in the previous section, the chemical potential in eq.( 1.50| ) is added to the Lagrangian as 



a Lagrange multiplier of a conserved charge. In this thesis we will consider the isospin chemical potential 
related to the isospin number, which indicates the difference of the baryon number between the u and the d 
fields. The isospin density number operator is defined as 

ni=\(y}u^Sd)^q^''-^q (1.60) 
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where q = (u, d) in the case of two flavors. Then the QCD-Lagrangian with massive quarks acquires the form 



2 

Cqcd{M,Q,\)x,t3U,Q). (1.61) 
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Chapter 2 



Tree-level masses, thermal 
propagators and vertices at finite 
isospin chemical potential. 



For radiative calculations, we need the different propagators derived from the free Lagrangian. To obtain 
the effective Lagrangian at the one-loop level, first we have to calculate the vacuum expectation value of the 
field U, named U. Taking the 0(p2) La grangian and setting tt" = we have that the effective potential is 

Veff = -^Tr [1^3, tj] [1^3, i/t] + 2BM {U + U^)] . (2.1) 
Taking the general expression for U 

U = cos V3 + «T • n sin (/J, (2-2) 

with 

n = cos0sin^^el + sin sin 6*62 + cosflea, (2-3) 

the previous expression becomes 

Veff = -f (i^f sin^ V3sin6' - m^cosv?) . (2.4) 
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Minimizing the efTcctivc potential, we find that the vacuum expectation value of the U fields depends on /ii 



[KST+00, KTOl and is given by 



U : 



+ i Ti cos ( 



T2 sm ( 



1 - 



if fif < w? 
if > 



(2.5) 



Now, we can expand the chiral Lagrangian in terms of the pion degrees of freedom up to the necessary 
powers for one-loop radiative corrections. The two phases, < m and |jUi| > m, will be referred to as the 
first phase and the second phase, respectively. The expansion of the chiral Lagrangian is explained in App. 



A.l. 



2.1 First phase 



For the case of the first phase, the relevant Lagrangian for one-loop corrections is Ci-ioop = C2.2 +'C2,4 + jC4 



as was indicated in Eq. ( 1.34 ), where the d^j terms are 



^,2,2 

^2,4 

-^4,2 



(5^0) 



2 ™2 2 



1 

2 

1 m2 4 1 
1 



-4|9/7r|'|^|' + (a|^|')Vm2(|7r|')' 



m 



6/2 

2 r 



2 |9/7r|' Til - 2 {dn^Y \ti\' + dnl-d \tt\' + m^nl \ti\ 



2 , 22 I |2 



2/4 19/^1 - 2to'(Z3 + h) \A + h {dTTo) - m\h + h- etdhWo 



using the standard definition of the pion fields 



(2.6) 



(2.7) 



(2.8) 



7ro=7r3, 7r± = ^(tti =F Z7r2). 



(2.9) 



In the previous expressions 



K|2 = TT+TT-, l^zTTp = (9+71+) (c?," TT" ) , 



(2.10) 
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and the derivative di± that appears in the Lagrangian is defined as 



df = ifiiu. 



(2.11) 



This definition of the covariant derivative is natural, since we know [Act85| that the chemical potential is 
introduced as the zero component of an external "gauge" field. The term in Eq. |2.8| is the quark masses 
ratio (m„ - m^)/ (m„ + m^). 

To obtain the DJp, first consider the equations of motion in momentum space, extracted from the free 

/ .-1 . n \ 



Lagrangian £2,2 hi Eq. (2.6) 



The non-zero elements of the free propagator matrix are^ 



^00' 

























(2.12) 



iD ^(j))oo = p^-m^ 
iD'\p)_^ = iD-\-p)^^. 



(2.13) 
(2.14) 
(2.15) 



It is easy to derive the tree-level masses, because the propagator matrix has no crossed terms. Then, the 
equations of motion act independently on each pion field. The previous equations are the equations of 
motion in momentum space iD^^(p)ij respect to the field TTj. From these equations, it is easy to associate 
the different tree-level masses, obtaining 



TTiQ — m. 



m± = m =p /ij. 



(2.16) 



where the the different pion masses are defined as their rest energy. 

Since our calculation will be at the one-loop level, we do not need, as was said previously, the full formalism 
of thermo field dynamics including thermal ghosts. Then, by the general formula in Eq. ( 1.5tj| ), the DJp are 



^Note that D ^{p)ij is different from 1/D{p)ij in this case because D+- and -D-+ arc anti-diagonal. 



18 



Fi gurG 2.1: Propagators in the first phase (and in second phase for ^ m). The double line denotes the neutral propagator. 
The single line denotes the charged propagator with an arrow pointing from + to — . 



Ddj{p) 00 

Ddj{p)-+ 



D^_{-p-T). 



(2.17) 

(2.18) 
(2.19) 



Diagrammatically, the neutral pion propagator will be denoted as a double line. The ± propagator will 



be denoted as a single line with an arrow pointing in direction from + to — , as is indicated in Fig. 2.1 



The vertices extracted from £2.4 and £4,2, shown in eqs. ( [2.7| ) and (2^), are 



yOO± 



^2A 



yOO 



3/ 



'—{2pq + 2p^p_- [p + q][p++p-] -|I^[pl-p°_] +?7i2-2//2| 



^{ - 2p+(7+ - 2p_q_ + [p+ + q+] [p_ + qJ\ 

+4/^1 [pI + q°-p°~ q°] + 2to' - S^i?} 

2 



2i- 



-{[- P+P- + -P-) + ^J'l]h - m'^ih 4-/4)}, 



(2.20) 
(2.21) 

(2.22) 
(2.23) 
(2.24) 



where p and q denote the external momentum of the tt" legs and p± and q± denote the external momenta of 
the TT* legs. In all cases, the momentum emerges from the vertex. 

Diagrammatically, using the convention for the propagators, a double line will denote a neutral pion leg 
and a single line a charged pions, where an arrow pointing away from the vertex denotes a 7r+ leg, and an 



arrow pointing into the vertex denotes a tt , as can be seen in Fig. 2.2 
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T/0000 T/°o± V±± 

^2,4 ^2,4 ^2,4 

[4^= — a — ^— 

^4,2 ^4,2 

Figure 2.2: Vertices in the first phase (and in second phase for |^i| > m). The double fine denotes the neutral pion leg. The 
single line represents the charged pion leg with an arrow pointing outside from the vertex for the 7r+ and into the vertex for the 
TT^ (or 7r+ and tt^ in second phase for |mi| ^ m) 



2.2 Second phase 

Let us denote the vacuum expectation value of the U field as 

U = c + ifis, (2.25) 

with 

-S' ..yTf, (...) 

where, for any vector, the v^ components refer to the components of a rotated vector by the azimuthal angle 

vi = vi COS0 + V2 sin(/), V2 = — vi sin + V2 cos 0, V3 = V3. (2.27) 
The vacuum expectation value of the U fields breaks SU(2)-^U(1). 
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The expanded Lagrangian, keeping all the details, is given by 



•^2,3 

^4,1 
^4,2 



^{[-(a7r)2+;,?(.2^2 + ^2_3,2^2) 

+2|^i|c(ifi9o^2 - ^290^0] TT^ + (tt • dnf^ 
/ifcs[4s2(/i + /a) - 2c2/3 - 5^/4]^^, 

-2s2 [8|Ai.|CTi<9o^2 + -3c2^2 _^2)](;^ 

+M?c2[s2^2 _ c27r2]Z3 + [c2(97r)2 + 2|//.|c(l - 3c2)iriao^ir, 



6/ 
1 

7 
1 

7 



(2.28) 
(2.29) 

(2.30) 
(2.31) 



(2.32) 



-M?c^{(.' 



Here, = M?c for further considerations, and I will consider a negative chemical potential (jUj = as is 



the case in neutron stars, where a condensate of 7r_ quasi-particles might exist [YP04|. The same happens for 
the 7r-)_, for the opposite sign of /i^. Note that, in the second phase, there vertices with three legs will appear 
from (£2.3)7 and with one leg from (£4.1). The latter is responsible for the counter-terms of the tadpoles, 
but they will not be considered in the approximation we will use next. 

The inverse of the free propagator for pions in momentum space, i.e., the kernel of equations of motion 
extracted from £2,2 is given by the matrix 



iD 



( 



\ 


























-2i\^ij\cpo 





(2.33) 
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To extract the tree-level masses, we calculate the poles of the determinant of D at zero momentum as is 

\iD-\p)\ =[pl-fi^][pl-fx^{l + 3c')]pl = 0. (2.34) 



prescribed in | KT01 ] 



We can identify the masses with respect to the results in Eq. ( 2.16 ), which have to match for fii = m. Then, 
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Figure 2.3: Pion masses as a function of the isospin chemical potential at tree-level 

the tree-level masses for the second phase are 



mo = m^. = |//i|VT+3c2, to_ = 0. 



(2.35) 



In Figure 2_^, wc can see the behavior of the masses at tree- level in both phases as a function of a negative 



chemical potential (see also | SS01 ] for the treatment at tree-level) 



We can identify the tts = ttq field, as in KTOl ], through its mass at |/ii| = m and because it is diagonal 
in the propagator matrix. Therefore, there will be no difhculties in handling this propagator. 



DDjip)a 



— fif + ie 



+ 2TmB{pa)S{p -Mi) 



(2.36) 



Unfortunately, this is not the case for the charged pions, which are mixed in a non-trivial way.^ 
The free propagator for the charged pions in momentum space is 



D 



P 



-2i\^ij\cpo 



(2.37) 



with Dij = Dij^^. The Dolan-Jackiw propagators can be constructed with the general formula indicated in 



^Although we can set (f> = 0, 1 will continue using the notation with a tilde to recall that the n fields are not the physical 
fields 
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Eq. (|L53t). 

The main difficulty with this matrix propagator is that it is very cumbersome to use it in the different loop 
corrections. This fact motivates to proceed in a systematic way, through an expansion in a new appropriate 
smallness parameter, namely s when > m or c when ^ m. 

It is not difficult to realize that the different vertices that appear in our Lagrangian will correspond to 
different powers of c or s, depending on the case. Although it could appear as a trivial correction, we will 
keep only the zeroth order in all calculations. As we will see, this procedure is not trivial at all and provides 
us with interesting information about the behavior of the pion masses as function of temperature and isospin 
chemical potential. 

If we scale all the parameters with in all structures needed for one- loop corrections, we have that: 



Propagators -> D(p;T, ^j;m) = ^D(p;T,c) (2.38) 

Mr 



h constants k{A) = [h - X - \n{cli^)] (2.39) 



Vertices K,m(p; A^/i /, "i. A) = -p^^^w„,m(p; A, c) (2.40) 

Integrals ^ / -^A*--^ ^ f -^K^-" (2.41) 

Corrections ^ S(po; T, /^,; /, 771, A) = — ^a(po; T, c; A) (2.42) 

where from now on the bar on a parameter means that it is scaled with It is possible then to expand 
the propagator in powers of s" for|yU/| > m and c" when 3> m. Doing the same with the vertices, the 
resulting radiative corrections will be expressed in powers of s and c. The problem with divergences is then 
also fixed since they appear in powers of c and s. 

As we said before, we will keep only the n = terms in the previous expansions. This approximation is 
non-trivial, since it allows us to explore the behavior of the renormalized masses precisely in the vicinity of 
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the transition point and also in the region of high values of the chemical potential. Further corrections could 
also be calculated, taking into account higher order vertices and propagators. 

We will avoid the region where c ^ s (i.e. arround c = s = -^) and consider only the region where 
c — > or s ^ 0. This excluded area becomes smaller when n (the order of the expansion) starts to grow. 
By demanding that 0"+^, s'"'^^ < ^i^)"^^: we achieve this condition. Note that if n = 0, we would exclude 
the region of the chemical potential where y^m^ ^ yu? ^5 VSm^. We remark again that by going to higher 
orders in our expansion, the excluded region becomes smaller, so this is the best bound. 



2.2.1 Propagators and vertices at order s° 

As was said before, the neutral pion propagator is diagonal with respect to the charged ones. Its propagator 
will be the same in both limits at any order. 

The propagator for the charged pions at 0(s°) is given by 



D 



-2i\^j,j\po 



(2.43) 



It is more convenient to work with the combination of fields 



(2.44) 



where the tilde tilde convention was described in equation ( ^.27 ) We remark that the 7f± fields do not 
correspond to the physical charged pion fields but to a combination of them. This fact is a consequence of 
the non-trivial vacuum structure in this phase. This combination is also not trivial due to derivative terms 
and the inverse D'Alembcrtian operator. I would like to mention that, in the first phase {fii < m), the '!t± 
correspond effectively to the charged pions and the propagator becomes anti-diagonal. Then, the inverse 
propagator for charged pions are 



D-_l Dzl J 



\ P ~ 2|^,|cpo - 2^ 



(2.45) 
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and the free propagator matrix is 




(2.46) 



The Dolan-Jackiw propagators for charged pions, including isospin chemical potential (enough for one-loop 
calculations), are 

Ddj{p)+- = 2 o/i ^. +'i^nB{\po\W~2\^i,\p^) + \o{s^), (2.47) 
P - 2|/ij|po + le /if 

Dd.j{p)-^ = Dojhp)^^, (2.48) 

and the other propagators, D_^^ and D__ are of order s^, so we will neglect them. Diagrammatically, the 
propagators in this limit are the same as in the first phase (see Fig. p?l| ). Note that, in the case where 
— 771, the fields tt^ can be identified as the physical fields n^. 



The relevant vertices at order s° are 



V,"y- = ^{2pq + 2p^p_-[p + q\[p^+p_]+\^i,\[pl-p'i]-fi^} (2.50) 

^2!/ = - 2p-<Z- + [p+ + <Z+] [p- + q-] 

-M^i^[pl +ql- / - q"] - 6p^} (2.51) 

2 

^4,2 = -2*y|{2p<7?4 + /i?(Z3 + ?4-ew'M} (2.52) 



4{[P^P- + \M -p")]h + ^^%}, (2.53) 



which, diagrammatically, are the same as the first phase (Fig. 2.2), where p and q denote the external 
momenta of tt^ legs and p± and q± denote the external momenta of the tt* leg. In all cases, the momentum 
emerges from the vertex. These vertices include higher terms of order s and s^, and other vertices also appear: 
V22, ''^2''4*'''i ''^J^** of ^(^^) ^^'^ "^2 3 ^ ^''-'Ts J ^4 1 of ^(*)- As was said, however, since we are working at 
zeroth order, these terms will not be considered in the analysis. 
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Ai 

1^22 

Figure 2.4: Propagators in the second phase for \fii \ 2> m. 



2.2.2 Propagators and vertices at order c' 



Proceeding in the same way as in the > m case, the propagator of the charged pions tTi and tTj at zero 
temperature at 0{c'^) in the region j^J ^ to is 



D 







(2.54) 



A was said before, the Doo propagator in Eq. ( 2.36| ) remains the same in all cases. The Dolan-Jackiw 
propagators for charged pions are 

Ddj{p)ii = DDjipU + ^Oic"), (2.55) 



1 



Ddj{jpU = ^— +27rnB(|p"|)5(/) + — 0(c2), (2.56) 



and the propagators and D21 are of order c. In the case of a propagator with zero mass, it is necessary 
to introduce a small fictitious mass as a regulator, i.e. D^^ = lim,,^o(p^ ~ ff )- Note that, in the chiral limit 
(c = 0), the fields tt^ and tTo have the same behavior. The other components of the propagator matrix are of 
order c. Diagrammatically, we will denote the D^^ propagator with a line and D22 with a dashed line, as we 



can see in Figure \2.% Note that when 00, the fields tt and tt can be identified as the physical fields 



7r+ and tt , respectively. 
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T/112 




4,2 



T/OOOO 

1^2,4 




Figure 2.5: Vertices in the second phase for 2> m). The double hnc denotes the vr'^ leg. The single line represents the tt^ 
and the dashed line represent the tt"^ leg. 



The relevant vertices at 0{c^) arc indicated in Figure 2.5 with the analytical expressions 



T/002 

•^2,3 



T/OOOO 

'^2,4 



T/0022 
''2,4 



T/1122 

'^2,4 



T/112 _ ^l/^iL o 
''2,3 — 'f~f(2) 



VI 



i 



3J2 + 2^(1)9(1) - [P(o) + 9(0)] [P(i) + 9(1)] + 4^^?} 

^{[2^(0)9(0) + 2^(2)9(2) - [P(o) + 9(0)] [P(2) + 9(2)] + 2ptj2| 
^{[2^(1)9(1) + 2p(2)9(2) - [P(i) + 9(1)] [P(2) + 9(2)] + 2Aif I 
^4,2 - -4*yi{p9^i +/9°^2 + + ^2)} 
-z4^[M + 2pV](^i+y. 



(2.57) 
(2.58) 
(2.59) 
(2.60) 
(2.61) 
(2.62) 
(2.63) 



These vertices include higher corrections of order c and c^. As in the > m case, other vertices of higher 
order in powers of c also appear: ^2^2: ^2%^ "2^3 j ''^2^3: ^2°4^i "^2^4^! "2^Fi ''^4^2 j of ^(c)- 
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Chapter 3 



Thermal pions in the first phase 

This chapter wih show how to compute radiative corrections to the masses, condensates and to the PCAC 
relation which ahows us to extract the pion decay constant. 

3.1 Masses 

In the previous section, we extracted from the effective Lagrangian the different vertices and propagators 
involved in the self-energy corrections. Using the fact that the anti-diagonal terms of the propagator matrix 
are the same, except that they have opposite flow of momentum, 

E(p)., = (3.1) 

we just need to calculate Eoo and S_+. The loop contributions to these self-energies are shown in Fig. 

The calculation of the different loop corrections is standard. For denominators of the form {k + buY — 
(or delta functions in the thermal insertions), we just change the variables of integration 4- hu^ k'^. The 




Figure 3.1: Loop contributions to a. Eoo and b. S^-i- 
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calculation of the different kinds of integrals can be found in appendices B.4 and 
Defining A = A — ln(TO^/A^), the resulting self-energy corrections are 



(3.2) 



)m.J 



= a{(p + /iiu)'[|A-2r4 + |(/o + /)] +8(po + /ii)e(Mi)? 

-m2[|A - 2h + \h - K^o - 2/)] }, (3.3) 
S(p)^_ = S(-p)-+, (3.4) 

where a = (m/47r/)^ is the perturbative term that fixes the scale of energies in the theory (for energies below 
47r/). The functions /, J and In are defined as follows: 



dx\/ — l[nB{mx — \fii\) + nsimx + |/ii|)] (3.5) 



/>oo 

J ~ dxx\/ x^ — \\nBijnx — — nsimx + (3.6) 



/>oo 

/„ = / - 1 2^5 (mx). (3.7) 



In all these radiative corrections, since the loops do not carry external momentum, it is not necessary to 
set p = 0. The resulting self-energy corrections are proportional to p^, po and constants. Proceeding with 



the renormalization of the corrected propagators, described in section 1.1.3, we choose the rcnormalization 



constants in such a way that they absorb corrections proportional to p^. Taking 



Zo = l + a[IX- 214 + ^1], (3.8) 
Z± = l + a[|A-2r4 + |(/o+/)], (3.9) 



we have that 



^r{p)oo = am^[- ^k~327ThlJj-2Io + il], (3.10) 
^r{p)-+ = a{(po + Mi)e(/^i)8mJ + m2[- i[3 + 2/o]}, (3.11) 
Sfl(p)^_ = S]fl(-p)_+. (3.12) 
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It is important to remark that radiative correetions will leave a dependence on the chemical potential for 
the pion mass only for finite values of temperature. In a strict sense, this procedure does not allow us to say 
anything new about an eventual chemical potential dependence of the masses at T = (cold matter) , which 
is already included in Ci. 

There are no significant difficulties in extracting the masses from the renormalized propagator matrix, 
since there are no mixed terms, so the determinant of the renormalized propagator is 



\Dr{v)\ = ^^'(p)oo D-^\v)-^ Dr{v)^- 



(3.13) 



For finite T and /x/, we find the following expressions for the masses : 



m^o(T,Aii) = m{l + a[-r3/4- 16^V,rfZ7-Io + 2/]}, 
m^+{T,fii) = m|l + Q;[-r3/4 + /o + 4e(^i)j]| -/ii, 
m^-{T,fi,) = m{l + a[-r3/4 + /o-4e(M:)j]}+/i:. 



(3.14) 
(3.15) 
(3.16) 



The term WTr^eud^h that appear in the m^o corrections, which value runs between 0.03 and 0.35 can be 
neglected. 

We can see that, at zero temperature and chemical potential, we obtain the usual pion mass 



= m[l - ah/A] w 140 MeV. 



(3.17) 



The different values for the tree-level constants can be found in App. B.2 



If the chemical potential of the charged pions vanishes, i.e for symmetric matter, at finite T we get the 
well known result 

m4T)=m{l + a[-h/4 + Io]] (3.18) 



of chiral perturbation theory [GL87a , see also |Lar86, CL90|. However, due to radiative corrections to the 
neutral pion propagator, its mass will acquire a non-trivial chemical potential dependence for finite values of 
temperature. 

Figure 3.2 shows the dependence of the tt" mass as a function of temperature and isospin chemical 
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0.5 



Figure 3.2: m^o as a function of T and fii in units of m,r- 



potential. It is a growing function of both variables, temperature and \fj,i\ 



From Figure 3.3 we see that, at zero temperature, we agree with the usual prediction, m^- = + fii. 
In fact, at zero temperature, the tt^ should condense when /ij = — (the inverse situation occurs for 
7r+). Now, this situation changes if temperature starts to grow. The condensation point no longer occurs at 
fii — —rriTr, but it moves to the left, i.e. to a bigger (in absolute value) chemical potential. In the vicinity of 
/.Ji = nij^ the mass starts to decrease. This figure shows the 7r~ mass in the normal phase only, < m^. 
Even though we can consider higher values for the isospin chemical potential in the first phase, since m > 



(see appendix B.2), this situation will be discussed in the next chapter. 

For the case of the 7r+ mass, the behavior is the same, exchanging for — /Zi. 



3.2 Decay constant 

In connection with the behavior of f-^[T^^,i) when /ii < TOtt, we have made use of PCAC, which provides us 
with a relation between the renormalized propagator and the pion decay constant. 

The axial current is obtained as the functional derivative of the action with respect to a"^, with = a°r°/2 

A;i = ^(A/,0,/i:4r3,0) (3.19) 



31 




Figure 3.4: Diagramatic contribution to the PCAC relation at one-loop corrections. 



The expansion of the axial current in terms of powers of the pion fields can be found in detail in Appendix 



A. 2.1. Defining = T the axial current is 



= -fdTr" (3.20) 

4,1) = -fi9i^)^ 

4i,3) - ^[2|7r|W-7r09|7rp] (3.21) 

4,3) = ^[2(^o' + 2Kn(9,7r)±-7r±a(7r2 + 2|7rn] (3.22) 

A^'gi) = -j2m^hdn" (3.23) 

= -j2m^h{dnr)^ (3.24) 



Saturating the axial current with a single pion (the reduction formula can be found in Appendix B.7), we 



32 



obtain 



ip^,f[^ + a[-jX + 2li 



]}^ 



i(p±^i'u)/|^/Z^ 
Tiu^fe{fii)a&mJ. 



a[ - |A + 2/4 



\iIo + I)]} 



(3.25) 



(3.26) 



The different contributions to the PCAC relation in the expansion of the axial eurrent are shown in Figure 



3.4. Using the values of the renormalization constants Za in equations (3.8) and ( |3.9D , the renormalized PCAC 
relation is 



{{0\A°\7r%{p))) = tpf[l + a[h-4l]}, 

{{0\A^\Tr]^ip))) = iip±iJ,u){l + a[h-2iI + Io)]}Tiueip,)8af 
Then, we define the effective decay constant as the part proportional to p^, so 



mj. 



/,o(r,/.,) = /{l + a[r4-4/]}, 
h±[T,y^i) = /{l + a[f4-2(/ + /o)]}. 



(3.27) 
(3.28) 



(3.29) 



We can see that, at zero temperature and chemical potential, we obtain the usual pion decay constant 



fn = f[l + al4] «92.5MeV. 



(3.30) 



For an increasing chemical potential, the couplings f.,^{T,pi) decrease faster (Fig. 3^). This effect is 



enhanced for /ttq (T, /ii) (Fig. 3_^.a) and is related to the fact that J^o (T, /ii) only receives radiative corrections 
from charged pion loops; in contrast, (T, pi) receives both loop contributions, from charged pions and for 
neutral pions. For the case of zero chemical potential, /7r(T) is again the same one for the three pions. 



/,(T) = /[l + a(r4-4/o 



(3.31) 
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3.3 Condensates. 



There are three kinds of condensates relevant for the present analysis: the chiral condensate {{qqj), responsible 
for the spontaneous chiral symmetry breaking; the pion condensate ((tt)), which gives us information about 
the pions in the condensed phase; and the isospin number density ((«/)), which tells us about the difference 
of u and d baryon number inside the thermal bath. The expansions of the currents are explained in detail in 



Appendix A. 2. In the case of the first phase, as is expected, the pion condensate is equal to zero. 



3.3.1 Chiral condensate. 

The chiral condensate can be defined as the expectation value of the scalar current in the vacuum, {{qq)) = 
((0| Js 1 0)), with the scalar current for massive quarks at finite isospin chemical potential defined as 

J. = ^(M,0,i/iiWT3,0) (3.32) 
The components needed for radiative corrections are 

((^.(1,0))) = -2Bf (3.33) 

((J.(i,2))) = B{{0\nl+2\nf\0)) (3.34) 

((^.(3,o)» = -^Bm^h + k + hi). (3.35) 

Using the fact that 

{{0\TT,{x)7r,{x)\0)) = ((0|r7r,(x)7r,(y)|0)),y^,, (3.36) 
the corrected chiral condensate turns out to be 

{{qq)){T, fi,) = -2Bf {l - a [\h - 2h ~ 32n^hi + 2/o + 4/] } , (3.37) 
where the chiral condensate at zero temperature and chemical potential is 



(qq) = ~2Bf [1 - a{^h - 2k - 32 A)] , 
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(3.38) 



and the chiral condensate at finite temperature and zero chemical potential is 

{{qq)){T) = -2Bf |l - a [ifg - 2k - 32TT^hi + 6/0] } 



(3.39) 



The constant hi, as was said before, is model dependent. If we accept the G-MOR relation 



(3.40) 



neglecting the term e^^^y that appears in the nij^o corrections, wc find from equations ( 3.1'i ) and (|3.30D that 



mlf^ = (m„ + md)Bf [l - a^^k - 2h) 



(3.41) 



Comparing with eq. (3.38), if we consider the G-MOR relation as valid, we can set hi = 0^. 



Figure 3.6 shows the chiral condensate at finite temperature and isospin chemical potential. It has a 



similar behavior as the pion decay constants in Figure 3.5 




Fi gurG 3.61 Quark condensate as a function of the temperature and the isospin chemical potential 



^For higher corrections to the G-MOR relation see |DFL96| 
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G-MOR relation at finite temperature and isospin chemical potential 



If we set hi = 0, ^ 0, from equations ( 3.18 ) and ( 3.31 ) we obtain 



m^{T)^U{Tf 



(3.42) 



Compared with equation ( 3.3g| ) we can see that the G-MOR relation is still valid. 

At finite chemical potential, however, this relation in not possible, since the states are not degenerate. 
This suggests to extend this relation using the average value of the results. 

First of all, we have to take into account that the tree-level part of this average cannot depend on the 



chemical potential. Observing equations (3.15) and (3.16), we can see that the only possibility is to use the 
average of the pion masses and the pion decay constants. Defining 



m^(r, /ii) 



i (T, pii) 4- m^- {T, pti) + m^o {T, pii)] 
m{l + a[-lh- ^32n\ljr + + 21)] } 



(3.43) 



/,(r,/xO = ^[f^^{T,^l,) + f,-iT,^l,) + f,oiT,^l,)] 

= /{l + a[f4-|(/o + 2/)]}. 



(3.44) 



and combining the quadratic terms, we find 



V(T, fi.ff^iT, ^i,f = (to„ + md)Bf {l + a + - ^327:hljj - 2/o - 4/] } 



(3.45) 



This result is exactly —^{rriu + ind){(qq)) in equation (3.37) but with hi = '^^uJt- We can say that the 



generalized G-MOR relation is 



(3.46) 



3.3.2 Isospin-number density. 

The isospin-number density condensate gives us information about the baryon number difference between u 
and d quarks in the thermal vacuum or populated vacuum. The isospin-number density is defined as the 
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0-Lorcntz and 3-isospin component of the vector current 



ni = \q^T^q = 



(3.47) 



We need to calculate the expectation value of the vector current in the vacuum. We will see that the only 
non-vanishing component of the vector current in the vacuum is the isospin-number density. 



Considering the expansion of the vector current explained in appendix A. 2. 3 and using the fact that 



kalo)) = 0, 



ka^3^3|0)) = 0, 



(3.48) 



we find that the only non-zero component of the vector current in the vacuum, for the first phase case, is 



:(V(i,2), 



i{TT+df_,TT_ - 7r_9^7r+) -I- 2/Xi|7rpu^|0))e3. 



(3.49) 



Using the fact that 



(3.50) 



(3.51) 



Figure 3/7 shows the isospin-number density expectation value in the thermal vacuum. It gets a non 
vanishing value only inside the thermal bath. 



If we compare the results of the PCAC relation for charged pions in equation ( 3.26 ), we find that the 
term proportional to is exactly the isospin-number density condensate 



(3.52) 
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Fi gure 3. 7; Isospin- number density condensate as a function of the isospin chemical potential for different values of the 
temperature. 
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Chapter 4 



Thermal pions in the second phase for 



> 



In chapter I presented a systematic expansion method for radiative calculations in the second phase. For 
the case of the chemical potential near the transition point between the two phases, it is possible to expand 
the radiative corrections in powers of s^: 

/, m; T, ^0 = «' ^ a„(p, T)s^^. (4.1) 

n=0 

If we consider the lowest order in this expansion, 0{s^), we find that the vertices and propagators are the 
same as in the case of the first phase. The only difference is that the factors m are changed by I remark 
that in this phase a negative chemical potential = — will be considered where the negative pion will 
condense. 



4.1 Masses 



Proceeding with the relevant vertices and propagators, the loop corrections to the pion propagator matrix 
are shown in Figures ^.l| and 2^ , respectively. 
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Defining A' = A — \n{fif /A'^), the self-energy is 



S(p)o 



a'{ - Ai? [|A + ifa - + 32^"^^ - utflj + 21'^ - |/'] 



-P [|A - 2h + |/' 



']}• 



«'{ - A^? [5^"3 + 327r2e2^;7 - 21^ - 8 J'] + [|A - 21, + + |/'] 

-2po|/^.| [|A - 2^4 + + |/' + 4J'] }, 



(4.2) 

(4.3) 
(4.4) 



These self-energies include also higher correction terms of 0{s'^). Now a' ~ (/ii/47r/)^ is the perturbative 
term, and the functions J', are defined as follows: 



/' = 
J' = 

/: = 



diV x'^ - 1 "^(iMiKa^ - 1)) + nB{\^J,l\{x + 1)) 



dxx V a;^ - 1 nB{\fJ.i\{x - 1)) - Jisd/iiKx + 1)) 



dx^/x^ - lx^"2nB{\ti,\x). 



(4.5) 
(4.6) 
(4.7) 



Note that these definitions are almost the same ones we used in the first phase. Here, however, the term 
that multiplies x in the argument of the Bosc-Einstcin distribution is fii instead of m. 



Following the prescription, indicated in section 1.1. 3| , that the renormalized self-energy does not depend 



on p , we have that the renormalization constants are 



Zo = 1 + a' [|A- 211 + §/'], 

Z± = i + a'[jX^2h + Plj + ll'], (4.8) 



and the renormalized self-energy corrections are then 



I]fl(po)-+ = a'{M?[-ir3 + 2/^ + 8J'] -8|m.|J'po}, 

Efl,(po)+- = Sfi(-Po)-+, (4.9) 
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plus higher corrections of order y.'fa' s^. We will use these results to extract the renormalized temperature- 
and chemical potential-dependent masses. 

To extract the to^o mass, we do not need any further effort, since it is already diagonal in the matrix 
propagator, and its renormalized self-energy is constant in the external momenta, i.e., it does not depend on 

m,o = + a'[ - \h - IGTT^lah -Io + 2/'] }. (4.10) 

For the case of m^± , the vanishing determinant of the charged matrix propagator, provides us with a 
second-order equation in powers of Pq: 



4 2 



2S]5,_^ + 4|m:|cS),_^ + + s'^i^^"j,^_^ + = 0. (4.11) 



The solution of this equation gives the tt^ masses, recognizing them respectively according to the tree-level 
case. We can write the expression for the masses in a more compact way, by expanding around the m+ term 
that appears in the last equation. The 7r+ mass is given by 

= m+ + aVil [ " ^3 + ^ - 4>/' + 0(s')] • (4.12) 

Some care has to be taken in the extraction of the 7r~ mass. For the case of m^o and m7r+ , since they 
have a tree-level mass, it is possible to consider them as 

mij(r, ^i) = rntini) [1 + a'a{T, s)] , (4.13) 

and, then, we can expand the corrections in terms of s (or c in the other limit). For the case of the tt^ mass, 
if we consider the expansion in mass corrections explained in chapter |l.l.3| , calculating the determinant of 
the propagator matrix, at the lowest order, we obtain a second-order equation for the Sm^ mass correction 
(which is the tt^ mass since it does not have a tree- level mass) . The result will be of the form 



Ai:|(5m_ = s'Ea + \ s^fi?^B + (s2Sa)2 + J^c^d + S|. (4.14) 
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Figure 4.1: m^o and m,^+ as a function of high vahics of the isospin chemical potential at different values of temperature. All 
parameters are scaled with rrin 



This mixes terms of order a' and y/ 0{a') + 0{a''^), so all these self-energy terms must be expanded inde- 
pendently in powers of s (or c). With this criterion in mind, the 7r~ mass is then 



(4.15) 



with 



S5t_+ = a'[-iZ3 + 24 + 8J'], 



(4.16) 



and Tc is the critical temperature, where the tt mass is zero at the transition point from the first to the 
second phase, i.e., when /ij = m: 



{Tc, to) = 0, Tc « 0.135m^ = 18.9MeV 



(4.17) 



To start the discussion, I would like to mention that, near the transition point, decreases, as in the 
tree-level approximation, and this behavior is enforced with temperature. In this region, to^o grows with 



both parameters (see Figure 4.1) 



For m^- , according to Figures |4.2| and we see that a certain critical temperature must be reached 
before the tt^ is removed from the condensed state. This is precisely the condition that determines the phase 



43 





.98 



1.02 



1. 04 



Figure 4.2: m^- (T, ^i). All the parameters are scaled with m-n- 




Figure 4.3: m^- plot versus isospin chemical potential for diflferent temperatures. All values are scaled with mT^. The vertical 
lino corresponds to Hi = m. 
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Figure 4.4: Phase diagram of the condensation point of temperature versus laoapiQ—chemical potential. The dashed hne 
corresponds to the approximation made by Sphttorff , Toublan and Veerbarschot 1 3TV02t ] . 



transition curve (Fig. 4.4). 

Starting from the first phase, according to the previous chapter, for T < myr, it is possible to find, at 
first order in T, an expression for the transition hne fJ-i{T) that coincides with the results given in eq. (8.91) 



from [3TV021:]. In Fig. 4.4, the dashed line corresponds to this approximation, valid at low temperature. 
However, this result gives us the transition line from the viewpoint of the first phase. 



The reader could think from Fig. 4.4 that the critical temperature remains constant in the second phase. 



However, the temperature actually rises as a function of the chemical potential, in the form 



transifJ-l) = — ~ 0.134|^i| 

m 



(4.18) 



The slope is very smooth, and therefore we cannot appreciate this behavior from the values shown in the 
figure. This growing behavior of the critical temperature is of course consistent with general statements 
about phase transitions in the Ginzburg-Landau theory and it has been actually calculated in the lattice for 
two or three color QCD. Nevertheless, if we consider the thermal corrections in the second phase, it happens 
that the condensation phenomenon starts to disappear for a certain value of T > 18.9MeV near the transition 
point, leaving a kind of superfluid state which includes the condensed as well as the normal phase (massive 
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pion modes). 



4.2 Condensates. 



As was said at the beginning of this chapter, the different propagators and vertices are almost the same as 
in the first phase if we consider the zeroth order in the radiative corrections (just changing m by |^i| and /ij 
by — /ii). In this case, the same occurs, as in the expansion of the different currents. 



4.2.1 Chiral condensate. 

Following the same procedure in the case of the first phase, the non-vanishing components of the chiral 
condensate at order s*' are 

((^,,(1,0))) = -2i?fc, (4.19) 

((J.(i,3))) = Bc40\TTl + m^\0)), (4.20) 

((J,(i,o))) = -ABfi^c{h + k + hi). (4.21) 

The term c must be kept because it is a global constant. 

The resulting chiral condensate at finite temperature and isospin chemical potential is then 

{{qq)) = -2Bfc{\ + a'[-\h + 2/4 - 2/^ - 4/'}. (4.22) 



In Figure 4.5, we can see the transition of the chiral condensate from the first phase to the second phase. 
Now, the chiral condensate in the second phase tends to decrease abruptly with the chemical potential. This 
effect is enhanced by temperature.^. 



^Due to the arguments of the previous chapter, we can set the constant h\ ~ 0. 
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Figure 4.5: Chiral condensate as a function of the isospin chemical potential for different values of the temperature. The 
vertical line denotes the transition point betwreen the two phases. 



4.2.2 Isospin-number density. 



The non-vanishing components of the condensed vector current at order s° are 



((^))(i.o) = -f\fl^s'e■,u, 

((V))(i,2) = ((0| - iinidn2 - n2dn,) - 2\^l,\\nfu\0))e3. 



(4.23) 
(4.24) 



Like in the case fot the calculation of m^- , we have to keep the tree-level part proportional to 



The isospin-number density condensate at finite temperature and chemical potential is then 



(4.25) 



In Figure \l.6[ we can see the transition of the isospin number density from the first phase to the second 
phase. Like the chiral condensate, in the second phase ((n/)) tends to decrease abruptly with the chemical 
potential. Once again this effect is enhanced by temperature. 
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Figure 4.6: Isospin number density as a function of the isospin chemical potential for different values of the 
temperature. 



4.2.3 Pion condensate. 

Now, in the second phase, the pion condensate is finite and give us important information about the condensed 
phase. The pion condensate is defined as 

((tt")) ^ ((0|g-z75r"9|0)) = ((0|j;|0)), (4.26) 

and is, like the other condensates, a QCD vacuum-state operator with the same quantum numbers of the 



pion field. From Appendix A. 2. 4, the non-vanishing components of the pion condensates are 



Pp(i,o)» = -2Bfse, (4.27) 
((Jp(i,3))> = i?c((0|7r2 + 2|7rno))ei (4.28) 
Pp(i.o))) = -^B^l^c{k + h)el. (4.29) 

As we did for the chiral condensate, we must keep the term s, since it is a global factor. We can see that the 
pion condensate is oriented in the direction ei, i.e., ((tt)) — ((7ri))ei. 

Proceeding in the same way as we did with the other condensates, the pion condensate at finite temper- 
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Figure 4.7; Pion condensate as a function of the isospin chemical potential for different values of the temperature. The vertical 
line denotes the transition point between the two phases. 



ature and isospin chemical potential is 

((^^i)) = 2Bfs{l + a'hifa + 2U - 21'^ - 4/']}. (4.30) 

It is important to recall that the pion condensates, like the chiral condensate, give an information about the 
quarks condensed in the vacuum; in this case the pion condensate is a mixture of u'^^d and d'^^u.'^ This 
condensate breaks parity. 

Note that this result is basically the same as the chiral condensate (except the global factors s and — c) 
if we neglect the term hi. Then the chiral condensate and pion condensate satisfy the relation 



s{{qq)) + c({ni)) = 0. (4.31) 
We will see that this relation is valid in the limit when | 3> m. 



Figure l.l shows the behavior of the pion condensate at finite temperature and isospin chemical potential. 



If we compare with Figure we can see that the chiral condensate diminishes with the chemical potential 
in a similar rate as the pion condensate. This means that the quarks in the chiral condensate mix together to 
form mixed u-d pseudo-scalar Cooper-pairs state. Due to the thermal bath, both chiral and pion condensates 



Although related, do not confuse the pion condensate with the number density of condensed pions. 



49 



decrease, as expected, due to the increase of the tt mass due to the thermal bath. 
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Chapter 5 



Thermal pions in the second phase for 



^ m 



In this chapter, we wiU finish the analysis of the in-medium radiative correction to the masses and condensates 
in the region of high isospin chemical potential. As in the previous section, we will expand the different 
corrections to the propagators and the condensates, but this time in powers of c. Nevertheless, we have to 
include a logarithmic term in the expansion since loop corrections give rise to terms of the form A — ln(/ij^/A^), 
which cancel with the terms A — ln(m^/A^) coming from the li coupling constants. The radiative corrections 
can be expanded as 

I](p;r,Ai:;m,/) = a'Y, [a„(p,r) + fT;-(p) Inc] c". (5.1) 

5.1 Masses 

As we did for the ^ m case, we use the relevant vertices and propagators to compute the self-energy 
corrections. In our approximation (©(c*^)), the corrections to the Dga and are the same. This happens 
because the vertices and propagators used in the loop corrections differ in quantities of higher order than c. 
The loop corrections arc shown in figure ^.l[ .a, for D„o (which are the same as those of Dn, exchanging 



the double line for a single line), and figure 5.1.b for D 



Note that, in this region, tadpole diagrams appear, which are absent in the previous case. However, at 
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Fi gurc 5.1! Loop contributions to a. Sqo (or Sn) and b. S22 at order 



the order c'', it turns out that these tadpoles vanish, because the vertex is proportional to pa, and the tail 
of the tadpole does not carry momentum. At the order c, non-vanishing tadpole diagrams will appear (and 
also new £4 corrections). Nevertheless, for high values of the chemical potential, the leading behavior will 
be given by our approximation. 



The self-energy corrections are 



S(po) 



S(po). 



00 



11 



«Vf{(P^ - 1)[|A' + 21nc- ffi - Ik + l + ll'o + l^'T^] 
+ 1 - 4/^ + 2Aiipo) + 16So(po) + 16Si(po) - 327rzCi(po)}, 

S(po)oo, 



(5.2) 



(5.3) 





(5.4) 
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plus corrections of order a' jifc . The functions An, i?„, C„ are defined as 



^i(Po) 

^2(P0) 

BoiPo) 
BiiPo) 
B2{po) 
Ci(po) 



dx [SpgX^ + (pg — In [po(^^ — x) + x — i 



dxpQ In [po(2;^ — + 1 — ie] 
da;a;nB(|/ii|x) 



Po — 2po2^ — 1 + ie 



da; Vx^-lns ( I I a;) 



dx\J x^ — lnB{\fJ.i\x) 



Pq ~ 2pox + 1 + ie 
Pi 











Pl + l 




2po 




2po 




2po 


) 



C2(po) = [0{po - 2) + 0(-po - 2)]p2^(po/2)2 - lnBi\Po/2\)' 



(5.5) 

(5.6) 

(5.7) 

(5.8) 

(5.9) 

(5.10) 
(5.11) 



As we can see from the previous equations, it is highly non-trivial to identify the renormalized mass, as 
was the case in the region \ fii\ > m. Therefore, we need to expand the propagator, as was explained in section 



1.1.3, around the tree-level masses, identifying the term proportional to pg, or in this case {mt + 5mY . 

For the case of Sfloo, it is enough to evaluate Eoo with pg = m^o = mg + <5mg, because D~j^l^^ is diagonal 
with respect to the charged-pion propagator. We expand around mg and renormalize it, according to section 



1.1.3. Sflii and S]i?22 need to be evaluated for two values, 7n^± = rn± + Sm±. The expansion in mass 



corrections of the functions (5.5) - (5.11) is explained in Appendix B.6 



The renormalization constants for the different masses are (evaluated at these masses) 



Zo[mg] = l + 2a'|i[A + 21nc-;i-2Z2] -7 



Zi [m+] 
Zi [m_] 

Z2[m+] 
Z2[m-] 



38^2 



3 ^0 



(ImiI)}, 



Zo + Oia'c'), 



1 + lot 



+ 31nc-;i-2/2 + f + |7r^,. 



2„2 T-' 



16^4 T* 512^6 T'' 

5 A*F 21 ^lf 



34 p 



l + 2a'{|A[i-2r2 + 31nc+|-i^ 



6i^2l}, 

8A'io- 21^12}, 



1 + 2a'||Ari - 2/2 + 31nc + 2 - |/^, - SA'ggj, 



(5.12) 
(5.13) 

(5.14) 
(5.15) 
(5.16) 
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where the mass inside the brackets is the mass around which the self-energy expansion was computed. The 
function Knm is defined as 

r d. f^'^-^\^,t i. (5.17) 



The self-energy corrections, associated with the different tree-level pion masses, as functions of the corre- 
sponding renormalized masses, are 



Si?(m,o)oo = Aa'\^l,\[[6~^^r^^ + 13I'„ + 2K2l 

+167rz^nB(|/i,|)] (m.o - + 3/^}, (5.18) 

S_R,("^7r+)ll = i;fl(?71^+)oo + M?C'("'c^); (5.19) 

SK(m,-)ii = -2aV'{5 + W'^'5 + 16/o-8/;-2i^2i}, (5.20) 
Si?(m,+ )22 = 2aVi|{[-|-^ + 16/vio-5A'n]|/ii| 

+ [-§ + |^ + 20A'n + 3ifi2]m,+ }, (5.21) 

^Rim^-U = ^J.^Oia'c^). (5.22) 



Note that the self-energy corrections actually have the form S]ij(m^i) = I]5j[TOi] + 'E^j^[mi\mT^i presented in 
Section [l.l.3| . 



Since we have already expanded the self-energy corrections in powers of the mass corrections {6mi)", we 
can neglect higher terms in the determinant of the propagator, finding the solution for Snii from the pole 
condition. For m^o and , the corrections Sttiq and 5m_|- are of order fiia' (neglecting terms of order a'c). 
For m^- , because there is no finite tree-level mass, we need to keep terms of the order a'^ in the expansion; 
i.e. we need to consider, in principle, (Sm^)'^ ^ ^lu'^/"^ . Remembering that SJ? is of order /i^a', and SJj of 
order /^la', we can expand the propagator, neglecting higher order terms. 



The resulting expressions for the renormalized masses are surprisingly simple. Many terms vanish, in- 
cluding also some complex contributions. 
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2.4 




1.7 1.8 1.9 2 2.1 2.2 



Figure 5.2: m^o and r?i^+ as a function of high values of the isospin chemical potential at different values of temperature. 



m,o = \fi,\[l + 6a'Il,] (5.23) 
= \fi,\[VTT3^ + 6a I'„] (5.24) 
m^- = ^iC'(a'c^). (5.25) 



Note that rtij^- vanishes once again in this region, i.e., the pion condenses again, in spite of the fact that 
the thermal corrected mass started to grow near the phase transition point. This behavior, however, could 
be just a fictitious result from our expansion up to order c°. 

In this region, m^o increases monotonically both with temperature and chemical potential. mj^+ , as the 
temperature and the chemical potential rise, becomes asymptotically close to m^o, as was expected (see 



figure. 5^). In contrast with the first region, the tt^ mass grows with temperature, and a crossover occurs 
somewhere in the intermediate region of the chemical potential for the temperature dependence, since near 
the phase transition point we have (Ti, /ii) > (T2, /ij), and this behavior changes in the high chemical 
potential region in such a way that m^+{Ti, fij) < m^+(T2,^i). 

55 





0.25 


— 1 


-1 — I — I — 1 — I — I — I — 1 — I — I — 


— ■ — 1 — ■ — 




1 — ■ — ■ — 




T — ; 










T = 










0.225 


V 






T = . 


5inn- 




■ 




. 2 








T = . 














^ ^^^^ 














0.173 


































■v. 














0.15 


















0.125 






■v. "■^^.^ 












0.1 


1 
















2 


2.2 2.4 


2.6 


2 


.9 


3 



l^r I /mjT 

Figure 5.3: Chiral condensate as a function of the isospin chemical potential for different values of the temperature for high 
values of the chemical potential. 



5.2 Condensates 



5.2.1 Chiral condensate 



The non- vanishing components of the chiral condensate, neglecting higher corrections of order c^, are 



((^.))(i.o) 

((^s)>(l,2) 
{{Js)){3,0) 



-2Bf'c, 



(5.26) 
(5.27) 
(5.28) 



keeping the global term c. Following the same procedure as in the previous chapters, the chiral condensate 
at finite temperature and isospin chemical potential is 



m) = -2Bfc ll + a' 



+ 327r2/ii +21nc-4//, 



2 /ttT 

Mi 



(5.29) 



Figure |5.3| shows the chiral condensate behavior at finite temperature for high values of the chemical 
potential. It continues decreasing with both parameters. 
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Figure 5.4: Isospin-numbcr density as a function of the isospin chemical potential at different values of the temperature in the 
high chemical potential limit. 



5.2.2 Isospin- number density. 



The non- vanishing components of the vector current, within the thermal vacuum, neglecting higher corrections 
of order (? , are 



((^M»(i,o) = -l/iil/'s'i^Mea, (5.30) 

TO)(i,2) = \vi.\s'mi^l^~i^\mu^B^, (5.31) 

((^M»(3,o) = -4|/i,|3s2(;^ + Z2)y^e3. (5.32) 

The isospin-number density is then 

ini)) = -ImiI/'s' {1 + 2a' {\h + \h + 21nc- 4^] } . (5.33) 



Figure 5^ shows the isospin-numbcr density for different temperatures and high values of the isospin 
chemical potential. Due to the logarithmic term, it starts to grow again with the chemical potential. As 
in the case of the 7r+ mass in this high chemical potential limit, it grows with temperature, in contrast to 
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2.2 



2.4 2.6 



2. a 



Figure 5.5: Pion condensate as a function of the isospin chemical potential for different temperatures in the limit of high 
chemical potential. 

the ^ m case. A crossover of the different temperature lines must occur somewliere in the intermediate 
region of the chemical potential. 

5.2.3 Pion condensate. 

The non- vanishing components of the pseudo-scalar current, within the thermal vacuum and neglecting higher 
corrections of order , are 



(Jp(l,0))) 
(«^p(l,3))) 
(«^p(l,3)» 



= -Bs{{0\7r'o+nt+ni))e, 



(5.34) 
(5.35) 
(5.36) 



The pion condensate is then 



= ^Bfs { 1 



Z4 + 21nc-4/^- ( — 

Mi 



(5.37) 



Figure |5.5| shows the pion condensate at finite temperature and high values of the chemical potential. 
Like the isospin chemical potential, it has a deflection and starts to decrease, due to the logarithmic terms. 
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Comparing with equation (5.29), if wc neglect the term hi, the pion condensate and the chiral condensate 
foUow the relation s{{qq)) + c{{tti)) =- 0, as was the case in the region > to. 
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Chapter 6 

Conclusions 



In this thesis, different aspects of pion dynamics within a thermal bath at finite isospin chemical potential 
were developed in the frame of Chiral Perturbation Theory and Thcrmo-field Dynamics. The extraction of 
the pion masses, pion decay constants, and different condensates were analyzed under these conditions in 
order to get a physical insight about the pion properties in a hot and dense pion gas. 

First, an expansion criterion was developed to calculate radiative corrections when the chemical potential 
becomes bigger than the pion mass (second phase), separating the analysis for the two cases, \fj,i\ > m and 
|/ii| 3> m. Nevertheless, in our analysis we restrict ourselves to < to,, in order to avoid strangeness effects. 
We also considered the so called first phase |/ii| < to, finding also corrections to masses and decay constants. 

Then we calculated the three pion masses at finite temperature and isospin chemical potential. Each of 
the three masses behaves differently under the introduction of the isospin chemical potential. For values of 
the chemical potential near the tree-level pion mass, one of the charged pions will condense (in this case, 
the TT^ when fij becomes negative). The condensation produces a symmetry breaking in the flavor group 
SU(2)— >U(1), inducing a two-phase structure. However, the thermal bath shifts the condensation point with 
respect to the tree-level case. It was also possible to extract the decay constants (which now are different for 
the TT^ and the tt*), the chiral condensate and the isospin number density. In the case of the masses, the 
decay constants and the chiral condensate, it was possible to re define and generalize the G-MOR relation 
for finite isospin chemical potential. 

The fact that the charged pions acquire a shifted mass at zero temperature is a special property of 
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the pions. Except for the charged pion masses, the isospin chemical potential effects appear only for finite 
temperature; this is also known as the "isospin chemical potential silver blaze problem"^ |Coh04| which 
indicates that at zero temperature nothing happens, unless the isospin chemical potential becomes bigger 
than the pion mass. 

Then, we proceed with the consideration of pion masses and condensates for isospin chemical potential 
values in the vicinity of the transition point. It was possible to find the transition where the tt" mass vanishes. 
For a certain temperature Tc, which depends on /ij, a normal phase and a condensed phase, as in a superfluid, 
coexist. With increasing temperature, the thermal tt^ mass grows and thus the pion condensate and the 
chiral condensate start to decrease. As the temperature and chemical potential increase, the isospin number 
density (absolute value) grows more rapidly than in the first phase. 

Finally, in chapter^, the pion behavior for the high isospin chemical potential values was studied, obtaining 
the pion masses and condensates. According to our results, the tt^ normal phase dissappears once again, but 
this is an artifact of the approximation involved. 

A further analysis should include aditional corrections due to higher order diagrams in our expansion 
parameter. The extension of the present discussion to the strange meson sector could also be interesting to 
pursue. 



From Arthur Conan Doyle's story, where Sherlock Holmes knew the identity of the thief because a dog doing nothing. 



61 



Appendix A 



Expansion in Terms of Pion Fields. 



This appendix summarizes the main steps in the construction of the moment expansion of the chiral La- 
grangian and currents in both phases. Here we wiU use the vacuum expectation value in the second phase 

U — eJ''^'^'^ ~ cos 93 + ifi siniy9 = c + ifis, (A-.l) 
where for any vector, the Vi components refers to the components of the vector rotated in an azimutal angle 

vi = vi cos + V2 sin (/), V2 = — vi sin + V2 cos (/), V3 = 173. (A. 2) 

For the case of the first phase we only have to set = (i.e. c = 1, s = 0) and </> = in all the results of 
this appendix. 

The unitary U fields can be expressed as 

U = Uo + iT-U = Uo + if-U, (A.3) 

which satisfy the unitarity relation Uq+U'^ = 1. On the other side, the field U is defined asU = U^^^e^^U^^^, 
with 

^ = r-^ = T-j, and U'^' = e'^'^^/^ = coa^ + ifisin^. (A.4) 
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The exponential can be expanded in terms of the ^ fields. Using the fact that ~ we have that 



n=0 ^ ' n=0 ^ ' 



Now putting the previous equation between the vacuum field J7^^^ 



U={c + ^ns) ^ + ^(e + ^.^1 + (c - ^ j^^^ (A.6) 



then 



n=0 ^ ' n=0 ^ ' 

i' - --.E^ + ft + (-%«-.)i:(tl2lL. (A.a, 

n=0 ^ ' n=0 ^ 

As we need the terms only up to C'(^^), the previous expression reduces to 



[/o = c-sl^-\ce + \slie + (A.9) 

IJ - (■■^ + ca-^^$'-^c6$" + ^.$^)ei + (|2e2+6e3)(l-^$') (A.IO) 
+ terms of O(C^). 



Another important thing to keep in mind is the fact that the terms i^i/ and R^n, vanish for the values of 
the external fields v = ^jiiTsu and a — 0. 
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A.l Expansion of the Chiral Lagrangian 



The lagrangian with massive quarks at finite isospin chemical potential is 



/:2(M,0,i/i:T3U,0) 
£4(M,0,i/i:T3U,0) 



^Tr [{Df,U)^Di'U + 2BM{U^ + U)] 

\h {Tr[iD^U)^D'^U]f 

+ ll2Tr [{D^U)^D,U] Tr [{D^U)^D''U] 

+ j{l3 + l4){Tr [BM{U^ + U)])^ 

+ \hTr [{D^U)^D>'U] Tr [BM{U^ + U)] 

-\h [Tr [BA/([/t -[/)])' 

+constants 



(A.ll) 



(A.12) 



with the covariant derivative 



DU ^dU ~ i^i^u[\T^, U] = dUo + ir ■ {dU - ^i,uU x 63) 



(A.13) 



First of all, we need to compute the traces. It is enough to calculate the traces of the main structures of 
the Lagrangian: Tr[{Dt,Uy D^U], Tr[M{W + U)], and Tr[M{W + U)]^. 



Using the shape of the covariat derivative in Eq. ( |A.13D 



Tr[{D^U)^D,U] = Tr 



^df,Ua - IT ■ [d^JJ - ^ijujy X 63)! |ft,.c., ^ ^ i/| 
2[a^C/oa^C/o + {d^U - i-hUf,U X 63) • {fi u)] 



(A. 14) 



Now, using the fact that 



BM = S[i(m„ + m.d) + i(m„ - md)r3] 
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(A.15) 



the other terms we want to calculate are 



Tr [BM{U'^ + U)] = Tr [2BMUn] = 2m'^Uo 
Tr [BM{U^ -U)]^ = Tr [~2iBMT ■ Uf = - We^^J/l. 



(A.16) 
(A.17) 



A. 1.1 C2 up to order tt^ 

Now we are gomg to expand the the structures that appear in £2, Tr[{Df^Uy D^^U] and Tr[A-I{U^ + U)], 



neglecting higher terms of 0{^^). From Eq. ( A. 14 ) 



Tr [{D^UYD'^U] = 2[{dUof + {dUf + 2^l,{Ul^oU2 - U2doUi) + ^lUui + Ui)] . (A.18) 



Using the results of Eqs. (|A.9|) and (A. 10) 



{dUof + {diif 
UidoU2 - U2doUi 



III + f/| 



2csii+c^ei+^i--s^e 



-icsiie^Uc^ii+i'2-s'e)e- 



(A.19) 



(A.20) 



(A.21) 



Even though we neglect constants and total derivative terms in the lagrangian. we will keep them in the 
previous equation for further considerations in £4. 

For the case of Tr[M{U^ + U)], from Eq. ( A.16| ), we just need to replace the results of Eq. (|A.g| ). 



A. 1.2 £4 up to order vr^ 



In the next part, the structures of the £4 Lagrangian will be expanded in order of apparition from Eq. ( A.12| ), 



neglecting higher terms of 0{^^). In this case the resulting constant terms and total derivative terms will 
be neglected, since they do not contribute. It will be used also the fact that 69oC2 = ^£,2dQ£,i + a total 
derivative. 
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Using the results of Eqs. ( A.18|)-(A.21) up to order 



Tr [{D,,U)^D^^U\ 



+0(C^) + tot. deriv. + const. 



(A.22) 



Using Eq. (A. 14) and inserting Eqs. ([A. 9]) and ( A.IO) up to order 



Tr [{D^U)^D,U\ 



= 8m?s2[2m?c.s|i + (9oO' + (96 



-0{£,^) + tot. deriv. + eonst. 



(A.23) 



From Eq. ( A.lfc) and Eq. (A.9) up to order 



Tr [BMiU'f + U)] 



4m'' 



= -4771* 2cs6 + c^r - s^ef +C'(r) + const 



(A.24) 
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From Eqs. ( |A.18| )-( |AT211 ), Eq. ( [A.16D and Eq. ( |AJ| ) up to order 



Tr [{D^U)'^D>'U] Tr [BM{U'^ + U)] 

= Am\dif + 2^i,{s^oi2 + c(liao6 - Mi)} 

= -4m2[/i2.s(l - 3c2)ei + 2mi(1 - Sc^^i^o^ - c(9^)2 

+/ifc{|52|2 _ (1 _ 3s2)|2 + II}] (A.25) 
+ tot. deriv. + const. 



From Eqs. (A. 17) and (A. 10) up to order ^ 



Tr [BM{U^ - U)] = -iM^UI 



(A.26) 



A. 2 Expansion of currents 

In this section the currents that wiU be used in this thesis to obtain information on condensates and on the 
decay constant wiU be expand in terms of pion fields. 
The currents are sorted in powers of momentum 



J = J(i) + >/(3) + • • ■ (A.27) 

For the case of the PCAC relation, we need to expand the axial current up to order tt'^ in ^(ij and order 
TT in A(3). As this relation will be used only in the first phase, it will be calculated only with the values c = 1 
and s = 0. So, the effective axial current will be of the form 



^eff = Aia) + Ai,3) +^3,1) (A-28) 

For the case of currents that will be used in the calculation of condensates, an expansion will be developed 
in the second phase as I did the previous section.. In the case of the first phase, we just need to take the 
terms c = 1 and s = 0. 
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For one loop corrections to the condensates, we need to expand up to order tt^ in J(i) and order zero in 
J(3) according to the power counting in chapter |1.1.2| . The effective current in this case will be 

"^eff = -^(1,0) + -^(1,1) + '^(1,2) + ^^(3,0)1 (A. 29) 

and then, the respective condensate is 

((J)) = ((0|J(i,o)|0)) + ((0|J(i,2)|0)) + ((0|J(3,o)|0)). (A.30) 



A. 2.1 Axial current. 



As was said before, we need the axial current in the first phase, to get, via the PCAC relation, the pion decay 
constant. Since it will be saturated with a single pion, we need to expand it up to order (tt)^. The axial 
current for massive quarks at finite chemical potential is 

^^ = ^(M, 0,1^1^X3,0). (A.31) 

The expressions for ^(i) and A(3) with finite quark masses and finite chemical potential are 

^(1)^ = ^Tr[^{ir^C/t}i?"t/ +/I.C.] (A.32) 

^{3)t^ khTr[i{^T'',U''}D^'U +h.c.]Tr[{D^UyD°'U] 

+ ^l2Tr[i{^T'',U^D(3US'^ +h.c., f3]Tr[{D''UyDf'U] 

+\UTr [^{ir^ U^]D^'U + h.c] Tr [BM{U^ + U)] . (A.33) 

First we need to calculate the trace of the following term 

Tr[i{\T'',U^]DU] ^ Tr[{iT''UQ + U''){dUo + iT -[dU - ^l^uU y.e3\]] 

= 2 [UdUo - UodU + fi.uUoU x eg]" . (A.34) 
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Using the expansion of the fields in equations (A. 9) and ( A.10| ) up to order with c = 1 and s = we get 
Tr [i{^T, U^}DU] = -2 m - Mi< x 63] [l - fl^] - ||(^ • d^- (A.35) 



This expression starts with terms of order ^. In the case of A(^^-^ , we need just expressions up to order ^, but 



from equation (A. 14) we can see that the term Tr [(D^J/)^!?^?/] is of order Then we can neglect the 
terms proportional to li and I2, and using eq. ( A.16| ) we obtain 



1(3) 



-2^m^ [d^ - fiiU^ X 63] 



(A.36) 
(A.37) 



A. 2. 2 Scalar current. 

We need the scalar current to calculate radiative corrections to the quark condensate (or chiral condensate) 
defined as 

m) ^ {{0\qqm ^ {{0\JM (A.38) 



The scalar current for massive quarks and finite isospin chemical potential is 



(A.39) 



The expressions for ^5(1) and J sCS) with finite quark masses and finite chemical potential are 



s(3) 



Js{3) = -^{h + h)BTr[U^ + U]Tr[BM{U^ + U)] 
-\l4,BTr[{D^,Uy D^'U]Tr[U^ + U] 
-AhiBTr\BM] 



(A.40) 



(A.41) 



For Js(2), expanded up to oreder (tt)^, we use the results of equation ( A. 9 ) 



^.(1) = -2/'SC/o - -2fB (c - s|i - \ce) 
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(A.42) 



Js{3) at order (7r)° is just setting the fields U ^ U. Then using the equation ( A.1S| ), with U ~ c + ifis 



(A.43) 



A. 2. 3 Vector current. 



The vector current will be used to calculate the isospin number density condensate, where 

im)) ^ ((0|igV3<z|0)> ^ mvim- (A.44) 
The pseudo-scalar current for massive quarks at finite isospin chemical potential is 

n? = ^(^^'0'5Mi"r3,0) (A.45) 

The expressions for V^")^ and Vj-g^ are 

^(1), = ^^,.h^[ir^C/t]i?A<^/ +h.c] 
^(3)m = 5^1^^ [5'^"' ^1 ^^C^ + /i-c] Tr [(D„C/)t£)"C/] 

+\hTr[~i[\T'',U^]D^'U +h.c] Tr[BM{U^ + U)\ 

First we need to calculate the trace of the following term 

Tr[~i[\T'^,U^]DU] = Tr[-i{U y.TY{dUa+iT -[dU - ^l^uU yie:i\]] 

= 2[U xdU - fi.uU x{Ux 63)] ° (A.48) 
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(A.46) 



(A.47) 



Expanding up to order 



2 
+2 
+2 



ei 
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(A.49) 



63, 



we see that we only need to replace it to obtain V^")^ up to order For the case of we only need to 



set S, — 0. Using equation ( |A.16D and ( |A.14| ), we get 



(A.50) 



A. 2. 4 Pseudo-scalar current. 



The pseudo-scalar current will be used to calculate the pion condensate, where 



j«975T q\' 



(A.51) 



The pseudo-scalar current for massive quarks at finite isospin chemical potential is 



(A.52) 



The expressions for and ^"(3) are 



ja 

•Jp(i) 



ja 



\fBTr [ir"(t/^ - U)] 

\{h + h)BTr [iT-^iU^ - U)] Tr [BM{U^ + U)] 
+ lUBTr [ir'^iU^ - U)] Tr [{D^UyD''U] . 



(A.53) 



(A.54) 



Considering 



Tr [ir°(C/^ - U)] = 2C/", 



(A.55) 
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then the expansion of Jpi^i-^ up to order eorresponds just to re placethe expansion of U°' (eq. (A. 10) and 
'^p(3 0) replacing the previous equation, equation (A. 16) and (A. 14) with ^ = 



Jp{i) — 25/^ 



r4i - ye) ei + 662 + + Oie) 



J, 



p(3,0) 



2Bs [2m'^c{l3 + k) + tifs'^U] ei 



(A.56) 
(A.57) 
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Appendix B 



Some useful formulas, identities and 
constants 



B.l li constants 



The li coupling constants are defined as the value of the corresponding couplings at the scale A = m, 
subtracting the divergent part proportional to A: 



li 

327r2 



7 1 , 

k+ln—-X 



(B.l) 



with 



A 



4-fi 



ln47r + r'(l) + 1 



(B.2) 



for i = 1,...,6. The table BJ was obtained from | Sch02 ; it shows the different li{A) and Ij constants 
evaluated at the scale A = m, the corresponding 7, factors and indicates also from which processes they 
where obtained. The other constants, hi, that appear in the chiral £4 Lagrangian, are model dependent and 
do not involve divergences. 
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Value 



-2.23 ±3. 7 pLSj ' 
-1.7 ± 1.0 pCG94 | 
-1.5 |BCE+97| 
-1.8 TCGLOlf 
-0.4 ±0.6 ]CGL01 | 



Obtained from 



TTTT D-wavc scattering lengths ©(P*) 
TTTT and Ki4 

TTTT D-wave scattering lengths O(P^) 
TTTT scattering ©(P^) + Roy equations 
TTTT scattering 0{P^) + Roy equations 



TTTT D-wave scattering lengths 0{P^) 

TTTT and Ki4 
TTTT D-wavc scattering lengths 0{P^) 
TTTT scattering O(P^) -I- Roy equations 
TTTT scattering 0{P^) + Roy equations 



1/3 



6.03 ± 1. 3 |GL84| 
6.1 ± 0.5 [|BCG94 | 
4.5 |BCE+97[ 
5.4 [ ICGLOTf 
4.3 ± 0.1 jCGLOlU 



2.9 ± 2.4 [pLSl 
\h\ < 16 jCGLOl l 



2/3 



SU(3) mass formulae 
/v;4 decay 



-1/2 



4.3 ± 0.9 jlGLSl 
4.4 ±0.3 |BCT9a 
4.4 ±0.2 icGLOlj 



scalar form factor 0{P^) 
TTTT scattering 0{P^) + Roy equations 



13.9 ± 1.3 |GL84t 
13.0 ±0.9 |BCT98| 



o\p^) 



-1/6 



h 



16.5± 1.1 |GL84| 
16.0 ±0.5 ±0.7 jeCTO 



(r2), 0(P4) 
vector form factor O(P^) 



-1/3 



5 X 10^^ F.L81 



7T^ — ?/ mixing 



Table B.l: L constants values. 



B.2 /, B and tud- 



The one-loop corrections for the pion mass and the pion decay constant are 



= m[l - a/3/4] w 139.57 MeV 
= f[l + ah] w 92.42 MeV, 



(B.3) 
(B.4) 



with a = (m/ATrf)'^. From these definitions and with the value of I3 and Z4 it is easy to derive the value of 
the tree-level pion mass and decay constant. 



a « 0.017 
m « 141 MeV 
/ w 86 MeV. 



(B.5) 
(B.6) 
(B.7) 



From the definition of the tree-level mass, 



Birriu + rrid), using the current-quark masses values we 
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can extract the value of B 



ruu + md = 5 to 11 MeV =^ B w 1.8 to 4 GeV (B.8) 
From the value of the light quark mass ratios, we can evaluate the tud constant 

^ = 0.2 to 0.7 =^ end = """"""^ = -0.2 to - 0.7 (B.9) 

rrid ruu + rud 

B.3 Functions and definitions 

The functions involved in the radiative corrections in the first phase are 

a = {m/Anff (B.IO) 
A = A-ln(mVA2) (B.ll) 



1=1 dx\/ — l[ns(mx — \fii\) + nB{mx + |/ii|)] (B.12) 



oo 



J = dxx\/ — \\nBijnx — |/ii|) — nsimx + (B.13) 



oo 



/„ = / dxx^'^\Jx'^ - 1 luBimx) (B.14) 



1 

being a the perturbative parameter. These integrals do not depend on the chemical potential sign, and grow 
with both, temperature and chemical potential 

In the case of the second phase, we denote the same functions with a prime: a', J', J', 1'^ which arc the 
same functions, but with |yUi instead of m. 

a' - (A*i/47r/)2 ^B^^g) 
A' = A-ln(/i2/A2) (B.16) 



/' = y dx^Jx^^ - \\nB{\iiA{x - 1)) + 715(1^1 IG^ + 1))] (B.17) 

/OO 
d.Ta;Va;2-l[nB(|Mi|(a; - 1)) - nB{\y.,\{x + 1))] (B.18) 

/OO 
dxx^'^^Jx'^ - 1 27iB(|/ii|a;). (B.19) 
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These integrals are also growing functions of the temperature but decrease with the chemical potential. 
Other intermediate function that appears in the self-energy corrections for ^ m is 



mn 



[x2 - {n/2y 



(B.20) 



B.4 Dimensional Regularization 



In the case of a loop formed by one free-propagator and without thermal insertions 
integrals that appear are 



the different 





i 




A 


d'^k 


ik^ 


{2n)<i 


A 


d<^k 


ikt'k" 


{2'n:Y 


A 



A 



(47r)2 

A2 
(47r)2 

^A">A2 

4(47r)2 



A 

1^X2 

A 



(B.21) 
(B.22) 
(B.23) 



In the case of a loop formed with two free-propagators and without thermal insertions: 
combine the denominators of the propagators using Feynman's identity 



0. 



first we 



AB Jo [xA + (1 - x)S]2 ■ 



then, the different integrals that appear are 

d'^k 



^4-d 
^4-d 
^A-d 



(2TrY [fc2 - A]2 
d'^k ik^ 
{2nY [fc2 _ A]2 

d'^k ik^'k'' 



1 



(47r)2 
2A 
(4^ 



+ 1-A 



i-A 
2 



{2Tr)d [fc2 _ A]2 2(47r)2 



(B.24) 



(B.25) 
(B.26) 
(B.27) 



Integrals proportional to fc'' are equal to zero. The diagrams with thermal insertions do not need to be 
regularized. 



76 



First phase 


Second phase 




Fo(|/.,|,±|M:|) = ^I' 




Fi(|a.,I,±I^,I) = ±^.r 


Fo{m,0) ^ 






i^2™(0,0) ^C(2m + l)T2™+2 



Tabic B.2: F„ evaluated with the different values that appear in the self-energy corrections. 



B.5 Diagrams with one thermal insertions. 



In the case of a loop formed by one thermal insertion (and without free-propagator): ^-w^ , the general integral 
that appear is 



Fn{a,b) 



(47r)2 



fc>i3(|fco-&|)<5(fc'-a2) 



(Ixx'^^^/x^ - l[ns(|ax - 6|) + (-)"riB(|ax + 5|)] . 



(B.28) 



Table |B.2| shows the function Fn with the parameters a and h evaluated at values that appear in the self-energy 
corrections. i^2m+i(i,0) and integrals proportional to ki are equal to zero. 



B.6 Diagrams with two Dolan-Jackiw propagators. 

The apparition of vertex with three legs in the second phase, in the high isospin chemical potential limit, 
produces loops with two DJp. 
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= «V? [(Po + l)A' + |p2 + 2Ai(po) 
i ' = 16a' fi^ Bq{po) 



(B.29) 
(B.30) 
(B.31) 

(B.32) 



with 



AiiPo) 

^2(P0) 

BoiPo) 

-82(^0) 
CiiPo) 

C2{Po) 








2a'fi^[pl{l-X')+A2{Po)] 



8a'fi^B2iPo) 



dx[3pQX^ + {pq ~ l)xj In [^0(2;^ — x) + X — ie 



dxpQ In [pq{x — a;) + 1 — ie] 

„.2 



dxxnBilfJ'ilx) 
dx\/ x^ — Insd/iija;) 
dx\/ x^ — Insd/iija;) 



Pq — 2poa; — 1 + ie 



Pq — 2pox + 1 + ie 

Pi 

pI - 2pox + It 



Pl-^ 




Pl-l 




Pl + 1 




2po 




2po 




2po 


) 



(B.33) 
(B.34) 

(B.35) 



(B.36) 

(B.37) 

(B.38) 

(B.39) 

(B.40) 

(B.41) 
(B.42) 



[eiPo - 2) + e{-po - 2)]pl^{po/2r - lnB(|po/2|)^ 
where po = Po/a*i- 

To extract the masses, we need to expand the self energy in terms of the mass corrections around the 
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G„(l) 


g;(i) 


g;:(i) 


^i(l)--f 


A[il) = -f 


^'i'(l) = -f 


Ml) --2 + ^, 


^2(1) = -2 + :^ 




So(l) = 


W) = 4(f)' 










B2{1) = i/vio 




= -A'lO - fi^ii - iA-12 


Ci(l) = 


Ci(l) = 


Cr(l) = 2^nB(|Mi|) 


C2(l) = 


C^(1) = 


C^'(1) = 



Table B.3: Zero, first and second derivative of the functions A„, B„, Cn, evaluated with po = 1 



tree-level masses. Then, we need to know the different values of the functions described in the last equations 



up to the second derivative in the energy, as was described in chapter 1.1.3 



with G„ = An , Bn or G„ . and frit can be 



(B.43) 



mo = 1, m+ = -v/l + 3c^, m_ ~ 0. 



(B.44) 



As we are working in the expansion at order c^, then Gi!^\m^ ) = Gjf^(l) + 0(c2). So we need to evaluate 



these expressions at po = and 1. Tables |B.3| and |B.4| show the different functions and their derivatives 
evaluated in = 1 smd pQ = 0, respectively. 
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G„(0) 


g;(o) 


Gr.(0) 




AUO) = 


^i'(l) = -1 


A2{0) = 


^'2(0) = 


A'm - 




So(0) = 


R'VnA 4 ( ttT\'^ 128 / 7rT\^ 




B[{0)^Q 


Bi'(O) = 2/^ - 10/( +8/^ 


S2(0)=0 


B'M = 




Ci(0) = 


C((0)=0 


c;'(o) = 


C2(0) = 


C^(0) =0 


C^'(O) = 



Table B.4: Zero, first and second derivative of the functions An , Bn, Cn, evaluated with po = 0. 
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B.7 LSZ reduction formula for the PCAC relation. 

Due to the fact that in the first phase the physical pions are identified as tt^ , tt* and the equations of motion 
do not mix these fields, the generalization of the LSZ reduction formula is 

((0|7r'^^(a;i)7r"^(x2)---^°"(x„)k''(p))) 

y"d4xe*^"((0|I?^T[^°i(.Ti)7r'^^(a;2)---^""(x„y(x)]|0)), (B.45) 

with Ui, 6 = 0,+,— and where V is the Euler-Lagrange operator 

= -d^-n? (B.46) 
p± = -d^^2i^i,d^^^i\~m^. (B.47) 

Then in the PCAC relation, the reductions made when we saturate the axial current with a single pion are 

((0|7r''(x)k^(p))) = (,5bo<5ao + <56±,5aT)e-*^'^ (B.48) 

((0|7r"n2:i)7r"Ha;2)^"na;3)|7r''(p))) 

= (-Jw-Jaio +<56±<5a,T)e~*^"M(0|T^"=(^2)7r"-'(x3)|0)) 

+ ('5w'5a.o +'5f,±5a.T)e"*''"M(0|T^"^(xi)7r'^M2^3)|0)) (B.49) 
+ (<560'5a30 + 56±<5a3T)e"'''"M(0|T^'^^(2:i)7r''^(x2)|0)). 
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